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In The Name of Allah, The Most Beneficent, The Most 
Merciful 


So Exalted is Allah, the True Sovereign. Do not be hasty with 
the Quran before its inspiration is concluded to you, and say, 
“My Lord, increase me in knowledge” 


(Allah Almighty is Truthful) 
My dear colleagues 


Searching for the best is the path of human interest at all times, and 
in our endeavor to provide the best and keep pace with the great 
scientific development witnessed in our contemporary world, we 
present to you this humble work entitled 


(Neutrosophic linear models and algorithms to find their 
optimal solution) 


Science is the basis for managing life's affairs and human activities, 
and living without knowledge is a kind of wandering and a kind of 
loss. Using scientific methods helps us understand the foundations of 
choice, decision-making, and adopting the right solutions when there 
are many solutions and many options. 


Through this work, we present a study of linear models using the concepts 
of neutrosophic science, the science that was built on the basis that there 
is no absolute truth, there is no confirmed data, issues cannot be limited to 
right and wrong only. There is a third state between error and right, an 
indeterminate, undetermined, uncertain state. It is indeterminacy. 
Neutrosophic science gave each issue three dimensions, namely (T, I, FP), 
correctness in degrees, indeterminacy in degrees, and error in degrees. It 
was founded by the American philosopher and mathematician Florentin 
Smarandache, in 1995 and came as a generalization of fuzzy logic that 
was founded by the scientist Lotfi. A. Zadeh, in 1965. Neutrosophics is a 
new science. It studies the different spectrums that a person can imagine 
in a single issue, which gives a more accurate description of the data of 
the issue under study and thus accurate results that leave no room for 
coincidence that help in making decisions that suit all the circumstances 
experienced by the work environment of the system under study. In our 


quest to search for the best. In this book, we present a study of linear 
models and algorithms to find the optimal solution for them using the 
concepts of neuroscientific science. We know that the linear 
programming method is one of the important methods of operations 
research, the science that was the product of the great scientific 
development that our contemporary world is witnessing. The name 
operations research is given to the group of scientific methods used. In 
analyzing problems and searching for optimal solutions, it is a science 
whose applications have achieved widespread success in various fields of 
life. The characteristic that characterizes this science is the development 
of mathematical models, tools and techniques that have the ability to 
express the concepts of efficiency and scarcity in a well-defined 
mathematical model for a specific situation. It has the ability to use 
scientific methods to solve complex dilemmas in managing large systems 
in factories, institutions and companies and helps decision makers in them 
to make optimal scientific decisions for the workflow. These issues were 
addressed according to classical logic, but the ideal solution was a 
specific value appropriate to the circumstances in which the data was 
collected. It does not take into account the changes that may occur in the 
work environment. In order to obtain more accurate results and enjoy a 
margin of freedom, we present in this book a study of neutrosophic linear 
models and algorithms to find the optimal solution for them. What is 
meant by neutrosophic models are models in which the data are 
neutrosophic values, that is, variables such as in the objective function, 
which expresses profit if the model is a maximization model, and 
expresses a cost if the model is a minimization model, which in turn is 
affected by environmental conditions. 


We take it in the form Nc; = c; + €;, where €;; is the indeterminacy, and 
it takes one of the forms ¢; € [Aja Ajo | or & E {Aja Aj2} or otherwise, 
which is any neighborhood of the value c;; that we obtain while adding 
The data on the issue then becomes the cost (or profit) matrix 
Ng = [c; aE gl, and also the fixed values that represent the right side of 
the constraint swings, which express the available capabilities of capital, 
time, raw materials, etc., and they are also affected. In environmental 
conditions, we take it from the form Nb; = b; + 6;, where oF is the 
indeterminacy of the required quantities. It can take one of the forms 


6; © [Mis Hig] or 6; © {His, Hig}, and the same situation applies to 
examples of variables in constraints that express quantities. The raw 
materials consumed in the production process are taken from the form 
Najj =aij+yVij, where y;j is the indeterminacy of the quantities 
necessary for the raw material i to produce one unit of product j. It can 
take one of the forms y;; € ne Pija|, or Vij € {Gija, Pijo}, which helps 


us obtain more accurate results and gives companies a margin of freedom. 
This book includes eight chapters: 

Chapter I: Study of neutrosophic linear equations 

Chapter II: Neutrosophic Linear Models 


Chapter II: The graphical method for finding the optimal solution for 
neutrosophic linear models 


Chapter IV: The simplex direct neutrosophic algorithm for finding the 
optimal solution for linear models 


Chapter V: The modified simplex neutrosophic algorithm to find the 
optimal solution for linear models 


Chapter VI: The simplex algorithm with a synthetic basis to find the 
optimal solution for linear models 


Chapter VII: Neutrosophic Conjugate Linear Models and the Dual 
Algorithm 


Chapter VIII: Some applications to neutrosophic linear models 


I hope to God Almighty that this work will achieve the desired benefit 
from its preparation. 
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Chapter I 
Study of Neutrosophic linear equations 


introduction: 


Through our endeavor to provide everything new and keep 
pace with scientific development, and given the great 
importance of the linear programming method as one of the 
methods of operations research, we found it necessary to 
reformulate the systems of linear equations and some of the 
methods for solving them using the concepts of neutrosophic, 
because all the research and studies presented using the 
concepts of this science gave more accurate results. Compared 
to studies presented using classical logic. 


1.1. Systems of linear equations according to 
classical logic: 


The sentences of linear equations in which the number of 
equations equals m and the number of variables in them equals 
n are given according to classical logic in the following general 
form: 


44X41 + Q42X2 Se we AynXn = b, 


A24X4 + A22X2 Se AgnXn = bs 


AmiX1 + Am2X%2 Fo $F AmnXn = Dm 
It is written in the following matrix form: 
A.X =B 
Where: 
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Aq1 2 + Ain by Xy 
a1 a22 -+d2n b X2 

A= B=| 7 X= 
Aam1 = Am2 + Amn bin Xn 


Where aj and b; are real numbers for all values of i = 


1,2,...,m and j = 1,2,...,n 

We distinguished three cases of the systems of linear 
equations: 

The first case: 

The number of equations is equal to the number of variables, 
ie.,.m=n 

The second case: 


The number of equations is greater than the number of 
variables, 1.e.,m >n 


Third case: 


The number of equations is less than the number of variables, 
ie., M<n 

Below we will present the systems of linear equations using the 
concepts of neutrosophic science, where we will take the real 
numbers aj; and b; as neutrosophic numbers, that is, of the 
form Nb; and Naj;, indefinite values. Perfectly determined, 
they can be any neighborhood of the real numbers a;; and b;, 


written in one of the forms. next: 
Nai; = Ajj + ij and Nb, = b; + Lj where ij E [Arij A2ij| or 
Ey; € {Ari pri i} or otherwise, then the systems of neutrosophic 


linear equations is written in the following form: 


ee 
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1.2. Systems of neutrosophic linear equations in 
which the number of equations equals m and 
the number of variables in them equals n: 
General form: 
NQ44X1 + NQ42X2 +++ NQqyXy = Nbz 
NQz4X1 + NAz2X2 + +++ + NAgnXyn = Nb 


NamynX1 + NAmnX2 + °° + NAmnXm = Nbm 


In the following matrix form: 


NA.X = NB 
Where: 
Nay, Na 2..Nain Nb, Xy 
NA = Nay, NQg2 ...NQgn NB= Nb» X= X2 
Nami Nang = Ndinn NB mn Xp 


To find the general solution to the previous systems of 
equations, we examine them according to the three 
aforementioned cases: 


The first case: 
The number of equations is equal to the number of variables, 
Le.,.m =n. 
We write the systems of equations as follows: 
NQ44X1 + NQ42X2 +++ + NQyyXy = Nbx 
NQz4X1 + NAz2X2 + +++ + NAgnXn = Nbz 


NQn1X1 + NAy2X2 +++: + NAnynXn = Nby 
oes 
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In matrix form: 


NA.X = NB 
Where: 
Nay, Najy2..Nain Nb, x4 
NA = Nay, Nag ...NQ2n NB= Nb» y= X2 
Nant Nang Nan NDp er 


The proverbs matrix is a square matrix whose determinant is 
Ay= |NA| 


Here we distinguish two cases: 


1- Ay= 0. This case results in two cases: 
a. If Ay= 0 and Ay, #0 where Ay x; is the determinant 
J 


resulting from the determinant of the matrix of proverbs 
Ay after replacing the column containing the proverbs of 
the unknown x; with the column of constants (the values 
on the side The second of the equations) then the 
systems have no solution. 

b. If Ay= 0 and Ay “= 0, this means that the systems of 


equations are not linearly independent, meaning that 
some of them are linearly related to each other. To 
address this case, we delete one of the two linearly 
related equations, thus the number of equations 
decreases and becomes m'where m'’=m-—1_ and 
m' <n, which is identical to the second case that will be 
dealt with later. 

c. If Ay# 0, that is, the systems of equations are linearly 
independent and the systems have a single solution, 
which can be found in several ways. In this research, we 


7 eo 


Neutrosophic linear models and algorithms to find their optimal solution 


study Gaussian- Jordan method, which is the basis for 
the direct simplex algorithm that we use to obtain the 
optimal solution for linear models. 


1.3. Gaussian- Jordan method for solving systems of 
neutrosophic linear equations in which m = n: 


To clarify the mathematical basis of this method, we write the 
equations in the following matrix form: 


Nay, Nay2 NQin X41 Nb, 


Nap, Nap NQon X2 = Nb, (1) 


NQn1  NQnz-»NQnn 1 LXn Nb, 
Or in the following abbreviated form: 
NA.X =NB (2) 
Since Ay= |NA| # 0, this means that the matrix NA has an 
inverse that is NA~?. We multiply both sides of equation (2) by 
NA~tand we find: 
NA71.(NA.X) = NA71.NB 
Hence, we get: 


I.X = NB' 
Which is written in the ee detailed form: 
10 0...0 Nb, 
0 1 0O.. -O} Nb, 
Sats ol) 48) 
0 0 0. 4 Nb, 


This process is the basis of ae Gaussian- Jordan method for 
solving a system of linear equations. In order to convert Figure 
(1) to Figure (2), we follow the following steps: 

1- We express Figure (1) in the following table: 
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Table No. (1): Table of equations 


ariables 
x1 X2 eo cee Xn NB 
Equations 
1 Nay, | Nay “ah Nain Nb, 
2 Naz, | Naz re Nazn Nb, 
n Nanz | Nan uf Nann Nb, 


2- 


Ul 
Nai; 


NDB} = 


We convert the matrix NA to the unit matrix I by 
processing the lines of the table so that we make all 
non-diagonal elements in all its rows equal to zero 
and the diagonal elements equal to one. When we 
want to remove the variable x, from the equation f, 
we follow these steps: 

We divide all the elements of row t in which we want 
to make x, equal to one by Na;;, so x, becomes 
equal to one and the other expressions change. 

We set all elements of the column with x, (except 
row t) equal to zero. 

We calculate the rest of the elements of the new table 
from the following two relationships: 


N N Na, ; _ NajyNays — NajsNag; 
tN Nd 
Nb, Nb,Nats al Na;;Nb; (4) 
(vb, — Nay Ne) = MN Noah 
; S Nats Nats 


The element Na;, is called the pivot element 


From the previous processing we get the following table: 
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Table No. (2) Final solution table 


Variables 
Nx, | Nx2 Jou NXy NB' 
Equations 
1 1 0 “aes 0 Nb, 
y) 0 1 ee 0 Nb; 
n 0 0 ae 1 Nb, 


Through the table, the systems of linear equations are written 
in the following matrix form: 

I.NX = NB' 
ce Nb, 


x2] _ | Nb, a 


o> 
oO 


0 0 0..1) Ln} Lb, 


Nx, Nb; 
NX> = Nb, 
NXy Nby, 


= Nx, = Nb, ,Nx2 =Nb,,...,Nx, = Nb; 
The second case: 
The number of equations is greater than the number of 
variables, 1.e.,m >n- 
In this case, we form a new system from the set of equations in 
which the number of equations is equal to the number of 
variables by excluding a number of equations of m — n. Then 
we solve the new systems as we did in the first case and 
replace the resulting solution in the equations that were 
excluded to ensure that they are satisfied. 
The third case: The number of variables is greater than the 
number of equations, i.e.,m <n- 


_- 
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In this case, we are faced with a set of equations of the 
following form, 
Nay4X1 + NQy2X2 +2 + NQyy Xm $7 + NQiyXn = Nb, 
NQ24X1 + NAg2Xq + °° + N@amXm + °° + NQonXn = hie 


(5) 
eeese % NG hs ver Nae tte op ek Na = = - NBp, 
Which is written in the following matrix form: 

x4 
x2 
Nay, =Nay2 «. NQq yy we Nain Nb, 


na a| N@21  N@22 »N@am--+-N@in| yp — |NBa| y — ee. 


Naing Naina ow Nn o  NQinn Nby 
Xn 
In the following brief form: 
NAgmn) -X(nay = NBaniy (6) 
1- We partition the matrix NA( 7) into two matrices: 
a- A square matrix of rank (m.m) and we denote it as 
NCoumm): 
b- And a rectangular matrix of rank (m.n — m) and we 
denote it NDimn-m) 
2- We partition the column matrix X(,,1) into two matrices 
X" may and X" m1): 
Then the systems of equations (5) are written in the 
following matrix form: 
X' (ma) 


NC IND pei | ail est 
[ (m.m) (mn m)| , Sua 


| =NBaniy (7) 
NConm): X' m1) + NDonn-m)- x (n—m.1) = NB om.) 
From them we find that: 
NConm): X' (m1) = NB omy _ NDmn-m): x (n—m.1) (8) 
Assuming that |VC| # 0, we multiply both sides in relation (8) 
by NC~? and we find: 


_9- 
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NC1.NC.X' = NC71.(NB — ND.X’) 
I.X'=NC71.NB—NC™1.ND.X) (9) 
Assuming that NC~*.NB = NB' and NC~*.ND = NDGnn-m) We 
find that: 


10 0..0] [*%1 

0 1 0O...0] | *%2 

0 0 0..14 L%m 
Nb, Ndi, Ndjp.. NG yas —m) Xm+1 
Nb} | Ndj, = Ndby .. Ndi as ie ; (10) 
ND Ndi Nees Naincn- m) 


It is transformed into a set of linear equations as follows: 


Nx, = Nby — (NdyiX%m41 + Ni2Xmy2 ++ + Ndi (n—m)Xn) 
Nx2 = Nbz — (Ndo1Xm+1 + Ndg2Xmy2 + + NG) (n—m)Xn) 


Nin = ND SN Ga tmaa tN Geman eo ENG one) 
This means that we were able to calculate m in terms of 
(n—m)in terms of (1N—™), Xm41)Xm42) ++»Xn,we note that 
the values of the variables x,, X2,..., Xm, it relates to the values 
taken by the variables X44,Xm42)--»Xn, Or in other words, 
what we give to the variables X44,Xm42)---»Xn, and that for 
every proposition of values such as B41, Pm+2) ---» Pyfor these 
variables we get a set of values for the variables x1, X2, ...,Xm 
1S: 

Nx, = Nb, — (Ndi, Bm41 + Ndj2Bm42 + °° + Ndi cq_—m)Bn) 
Nxz = Nby — (Nao; Bm41 + Ndo2Bm+2 + °° + Ndgq__m)Bn) 


NXm = Nbm vs (NdiniBm+1 + Ndn2Bm+2 5 a Ndin(n—m)Pn) 
Thus, we obtain a solution that includes all the variables of 
sentence (5) 


-10 - 
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The solution is arranged as follows: 

(fi, Bo, sony Bm Bm+1) Bm+2) ee Bn) 
But since the variables X%44,%m42)+-,X,can take an infinite 
number of qualitative values (even if they are restricted by 
certain conditions), we obtain an infinite number of 
corresponding values for the variables x1, X2, ..., Xm. 
Therefore, if |NC| # 0, then the set of equations (5) has an 
infinite number of acceptable solutions of the form: 

(X4, X20 Xm Xm4iy Xm+20 > Xn) 
Thus, we obtain a solution that includes all variables of the 
sentence, which is the ordered solution: 

(fi, Bo, nay Baas Bm+1) Bm+2) ony Bn) 
1.4. Basic solutions of the neutrosophic linear 

equations: 

Since sentence (5) has an infinite number of acceptable 
solutions, we will try to limit ourselves to a limited number of 
them by setting the variables X44,Xm42) ---»Xn equal to zero. 
Then sentence (9) takes the following form: 


10 0.0] [%1] [Nb 
0 1 0..0] |*2}_ | Nb; 


Xm) LINDY, 


(11) 

0 0 0s 4 
From it we get: 

KS ND: 9X5 NDS ee ye SND 
So, the complete solution is: 

(Nb,, Nb35, ..., NDyy, 0,0, ... 0) 

We call this solution the basic solution because it is attributed 
to the rule with single normal vectors in the space R™ as 
follows: 


ee ie 
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1 0 0 
0 1 0 
ey = 0 C2 = 0 Cm = 0 
0 0 0 
0 0 m 


The set of vectors e€;,@2,..-,€m form a rule because they are 
linearly independent, and the vector NB’ can be expressed in 
terms of it using the factorials x1, X2, ..., Xm aS follows: 
NB’ = 1X1 + €2X_. ++: + emXm 
We call the variablesx,, x2, ...,X.,, basic variables and we call 
other variables X44, X%m42) ++» Xn 
free or non-basic variables because they take qualitative 
values. 
The process of choosing the variables x1, Xz, ...,Xm to be basic 
variables is a random process, as we can form other basic 
solutions, knowing that the possibilities available to obtain 
basic solutions are: 
ae n! 
ea m!(n—m)! 
It is a finite number of infinite acceptable solutions. 
Example 1: 
Find the joint solution of the following two linear equations: 
2X, + 7X2 + 3X3 + 2x, = [2,5] 
3x, + 9x, + 4x3 +x, = [3,7] 
X, + 5x2 + 3x3 +4x, = [4,8] 
In the set of equations, the number of variables is n = 4 and 
the number of equations is m = 3. Therefore, the number of 
basic variables is equal to 3 and the number of non-basic free 
variables is 7 —-m=1. The number of possible solutions is 
calculated from the relationship: 


[2's 
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a n! 
oe m!(n—m)! 
re: 

ce eke = 4 
4 31(4—3)f 
Write as follows: 

(X4,%2,%3,0) ,(%1,%X2,0,X4) , (%1, 0,3, %4), (0, X2, X3, X4) 
To obtain these solutions, we write the systems of equations 
in the following form: 

2X, + 7x. + 3x3 = [2,5| — 2x, 
3X, + 9x2 + 4x3 + x4 = [3,7] — xX 
X, + 5x, + 3x3 = [4,8] — 4x, 
The previous sentence is written in the following matrix form: 


2 7 3) [%4 [2,5]]} [2 
E 9 4 |: = |[3,7] -|3}-t (*) 
1 5 3] 1%3 [4,8] 4 
(iis ae x4 [2,5] 2 
=|; 9 4) X'= |X2] NB = }[3,7]} D=]1| X =bal 
L253 X3 [4,8] 4 
We calculate the determinant |C| . 
We find: 
2 3 
ICl=|3 9 44=-3 40 
15 3 
To find the solutions, we find the reciprocal of the matrix, 
2-7 3 
C= E 9 i 
f-c5 3 
We find: 
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a7 ; =f 
3 3 
ee ; —1 
3 3 
= as | 


We compensate in the relationship: 


NC-1.NC.X' = NC“1.(NB — ND.X’) 


We get: 
—7 —1 
ie, Se. 28 
5 —1].]3 9 4 
aaa cae 
—2 1 1 
—7 5 —1 
ee 3 | /1[2.5]] 72 
—-|5 =) [3,7]] —}1]. [x4] 
ES Ga cee 
3 3 [4,8]| 14 
—2 1 1 
1 0 OFfNX% 
0 1 O;|Nx2 
0 0 1ILNx3 
—7 5 —1 —7 5 —1 
2 3°] [25]} | 3 3B 2 
=| 5 —1].}[3,7]]-—| 5 —1]}.]1]. [x4] 
—-— -1 — = -1 +||4 
3 3 1 [148] 3 3 
—2 1 1 —2 1 1 


It is transformed into the following systems of equations: 


eee 
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~ fel 


Setting the free variable : equal to zero we get: 
On 
73 


Nx, = F | 
73 


1.€., 


=1 4 
X41 = lo, =| ,X2 = E 5| X32 = [3,5] 
Thus, we obtain the first neutrosophic basic solution, which is: 


(X1,X2,X3,0) = ([o, =|- 1,5].1351.0) 


In the same way we obtain other basic solutions. 

Dissolved basic solutions: 

The base solution is a degenerate and invalid solution if in the 
final result we obtain a value of zero for the variables that we 
chose as the base. 


1.5. Gaussian- Jordan method for solving a set of 
linear equations in which m < n: 
Based on the previous mathematical principles, the basic steps 
of the Gaussian- Jordan method are as follows: 
1- We write the systems of equations (5) in the following 
matrix form: 


1.X'+NC7!.D.X = NC7!.NB = NB' 


AR 
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1 PO tee ' 
[,NC-1.ND]. a = NB’ (12) 
Which is written in the following detailed form: 


10 0..0 Ndi, Ndjz..Ndicq— vel ff te 
Nbi, 


01 0..0 Nd}, Ndi « ' Nan m) (13) 


0 0 0. 1 Nd'y Ndlyg dinin m) 


The transition from Figure (5) to Figure (12) is done by 
following the same steps that we mentioned in the previous 
paragraph, but this method does not give us a basic solution 
unless we set all the free variables equal to zero. If we do that, 
we only get the first solution, and to get all the solutions, we do 
the following steps: 


a. We organize the following table: 


Table No. (3) The first table for the Gaussian- Jordan method 


ariables 
: Xy | Xq | wee | Xm | Xm. | Xmi2 | | Xn | NB 
Equations 

if 

1 Qi1 | G12 | + | Gaim | Gimta | A1m42 | | Gin | Ndi 

Ul 

2 Az1 | 422 | «+ | Gam | Gam+1 | 4am+2 | + | Gan | Nb2 

m a a a a a a Nb, 
m1 m2 ais mm MmmMm+1 mm-+2 ss mn m 


b. We find the identity matrix Ij, by processing the rows 
of the previous table in the same way that was explained 
in the previous paragraph. This is done by specifying the 
variables that will be entered into the base and let them 
be x1,%X2,.-,Xm. AS a result of this processing, we get 
the following table: 
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Table No. (4): Table of the first basic solution 


ariables , 

: Xy | Xz | 00 | Xm | Xms1 | Xms2 | -- Xn NB 
Equations 

1 1 | 0/...] O | Nady, | Ndjs |... | Ndinoim | Nb, 

2 0 1)... | 0 | Ndé, | Ndj, | .. | Ndon_m | Nb 

m 0 | 0]... ) 1 | Ndi | Ndb.| ..| Ndinn-m | Nbo 


c. Setting all the free variables in Table (4) equal to zero, 

we obtain the following first basic solution: 
(Nb,, Nb3, ..., Nbyy, 0,0, ... 0) 

d. To obtain a second basic solution, we replace one of the 
basic variables, say x,,, with one of the non-basic 
variables X41, by selecting the appropriate pivot 
element, and here it is Nq,,,. We work to delete x44 
from all equations except equation m .We set the 
coefficient of this variable in this equation equal to one. 
We perform the appropriate calculations through the two 
relations (4). We solve the following second basic 
solution: 

(Nb,, Nbb, ..., Nbj,_ 1,0, NDjy44, 0, «..,0) 
To obtain other basic solutions, we repeat what was stated in 
step (d). 


ee 
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1.6. Non-negative basic solutions of systems of 
neutrosophic linear equations: 

If all or some of the variables are conditioned to be non- 

negative, then some of the basic solutions are unacceptable 

because they violate the condition. In such a case, we have to 

look for positive basic solutions from among the basic 

solutions. 

Due to the difficulty of applying the method mentioned in the 

example, especially in the case that includes a large number of 

variables, the Gaussian- Jordan method was developed so that 

positive solutions are directly obtained. 

The new method was called the simplex method, which is 

carried out according to the following steps: 

1.7. The simplex method for finding non-negative 
basic solutions to a system of linear equations 

in which m < n: 

In the systems of equations (5): 

1- We make all elements of the constant’s column NB 
on the second side of the equations non-negative, by 
multiplying the equation whose second side is 
negative by (-1) 

We put the coefficients of the new systems in a table. 


i 


3- We form a rule consisting of m variables by choosing 
the variable that we want to enter into the rule, for 
example, x,, and then we calculate the index. 

Nbj|__ Nb; 

Najs = Nats 


We call the element Na,, the pivot element, we 


6 = Min| > 0; Na;; > 0,Nb,; > 0 


delete the variable x, from all equations according to 
the Gaussian- Jordan method, except for the equation 
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t, in which its coefficient is equal to one. We repeat 
the previous step until we form a base consisting of 
m variables. 

4- Setting the non-basic variables equal to zero we 

obtain the following non-negative basic solution: 
(Nbi, Nb5, ..., NDip—1,0, NDin 4.4, 9, «.,0) 

5- To obtain other new non-negative basic solutions, we 
choose one of the variables to be a basic variable, 
then we determine the pivot element and repeat the 
work we did for the previously mentioned variable 
x,;. We obtain a new non-negative basic solution, and 
sO we continue working until we obtain all non- 
negative basic solutions. 


We explain the above through the following example: 
Example 2: 


X, — 3x4 + 2x, = —[1,3] 

X_ + 2x, — 3x5 = [2,8] 
We multiply the first equation by (-1) until the condition 
Nb, > 0 is met, we obtain the following new systems: 


—xX, — 3x3 — 2x. = [1,3] 


Xo + 2X4 — 3x5 = [2,8] 

The stopping criterion is when we do not find a free column 
that we did not use for switching that contains a positive 
element (at least one), meaning that all the elements of the free 
columns that were not used during the swap are negative 
values. 

In the systems of equations, the number of variables is n = 5 
and the number of equations is 
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m = 2 Therefore, the number of basic variables is equal to 2 
and the number of non-basic free variables is n — m = 3. The 
number of possible solutions is calculated from the 
relationship: 
bee n! 
as m!(n—m)! 
ions 
5 5! 
o> Gao 
Write as follows: 

(x1, X2,0,0,0) , (x1, 0, x3, 0,0), (x3, 0,0, x4, 0), (x1, 0,0,0, x5), 
(0, x2, X3, 0,0)(0, x2, 0, x4, 0), (0, x2, 0,0, x), (0,0, x3, X4,0), 
(0,0, x3, 0, xs), (0,0,0, x4, xs) 

To obtain these solutions, we organize the following table: 
Table No. (5): The first table for the simplex method 


Variables 
X41 X2 X3 X4 Xs NB 
Equations 
1 —1 0 0 3 —2 [1,3] 
2 0 1 0 Z —3 [2,8] 


To find a basic solution to the set of equations, we choose a 
variable, for example x,, to be a basic variable, and to 
determine the appropriate anchor element, we calculate the 


index: 
Nb; [1,3] [2,8]) [1,3] 
6 = Min |——| = Min |——_ ,——-| = ——— 
in| | n| a 2 3 
That is, the fulcrum is a,, = 3. By performing the necessary 
calculations to delete the variable x, from the two equations, 
we obtain the following table: 


is 
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Table No. (6) The second table for the simplex method 


Variables 
x Xz | X3 X4 Xs NB’ 
Equations 
x4 il 0 0 1 _ 2 E q 
3 3 Be 
2 2 1 0 0 5 i 6| 
3 3 3° 


We choose another variable to be a basic variable. We note that 
the variable x, is ready to be a basic variable, and thus we get 
the following table: 

Table No. (7): Final solution table 


xy Xz | X3 x4 Xs NB’ 

CaN il 0 0 1 2 E 1 
3 3 3° 

Xp 2, ||| ab |) $6 0 5 é 7 
3 3 cu 


Thus, we obtain a base consisting of the variables x2 ,x,. We 
set the free variables equal to zero, and we obtain the following 
non-negative neutrosophic basic solution: 


(0.54)-0.55-1)-) 


To obtain other solutions, we repeat what we did to determine 
the previous solution. 


a2 s2 
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Conclusion: 

As a basis for neutrosophic linear programming, we presented 
in this research a study of the sets of neutrosophic linear 
equations, and the Gaussian- Jordan method, which is 
considered the mathematical basis for the simplex method used 
to find positive basic solutions, which can be used when there 
are restrictions on some or all of the variables to be positive 
values, which in turn was the basis for the method. Direct 
simplex used to find the optimal solution for linear models. 
Through the examples that we presented on the systems of 
neutrosophic equations, we arrived at basic neutrosophic 
solutions that express unspecified values. Such sentences can 
be used in cases where the data provided to the systems that 
operate according to these systems of equations is subject to 
change. Here we can benefit from the margin of freedom 
offered by neutrosophic values. 
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Neutrosophic Linear Models 
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2—1- Basic formulas of neutrosophic linear models. 


2-1-1-The general formula for the neutrosophic linear model. 
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model. 
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Conclusion. 


roe 


Neutrosophic linear models and algorithms to find their optimal solution 


Chapter II 
Neutrosophic Linear Models 


Introduction: 


In this chapter, we present the formulas of neutrosophic linear 
mathematical models, by which we mean linear models that 
contain in their mathematical relationships neutrosophic 
values, whether in the goal-function relationship or in the 
constraint relationships, and which take into account all the 
changes that may occur in the operating environment of the 
system represented by the model, which It ensures that the 
facility has a safe workflow, meaning that we will take the 
variables in the objective function as neutrosophic values, 1.e. 
Neo = ce; + & 

Also, the values that express the available capabilities are 
neutrosophic values, i.e., Nbj = bj + 6; and , Nay; = aij + Mi; 
where (Vj = 1,2,...,n,i = 1,2,...,m) are undefined values that 
have a margin of freedom and are taken according to the nature 
of the situation represented by the linear model, then We 
present the basic formulas of linear models through the 
following study: 


ee 
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2—1- Basic formulas of neutrosophic linear models: 


Neutrosophic linear models can be classified according to the 
following formulas: 


2-1-1-The general formula for the neutrosophic linear 
model: 


The general neutrosophic formula for the linear mathematical 
model is given in abbreviated form as follows: 


n 
NZ = XG BE Ej )x; — Max or Min 
j=l 


Constraints: 


n > 
>, Naix; (<) b +6; ; t=1,2,..,m 
j=l = 


Where cote, b +6, ay jf =1,2,..,n,i=1,2,..,m 
are constants having set or interval values according to the 
nature of the given problem, x; are decision variables. 

It is given in the following detailed form: 

Find 


NZ = Neyx1 + Nc2X2 +++: + NeyX, — Max or Min 
Constraints: 


= 
Nay X1 + NQj2X2 +++ + NAinXn (=) Nb; i=1,2,..,m 


Kip Xo, angry 0 


95%. 
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Linear models can also be expressed using matrices, and 
therefore the neutrosophic linear model given in the general 
form can be written using matrices as follows: 


Find 
NZ = NC X — (Max or Min) 


2 
NAX (<) NB 


Constraints: 


X20 
Where: 
Nay, Nay42..Nayn Nb, Nex x4 
NA= Nay, Naz ... NA2n NB = Nb, NC Ncp X= X2 
Nam Naqz NO Nb Newt — Lxn 


2-1-2- The canonical neutrosophic formula for the linear 
model: 


The linear program is set canonical if all the variables are 
constrained to be non-negative and if all the constraints are 
given in the form of inequalities must be placed in the form 


( <1is less than or equal to). The neutrosophic canonical form 
is written in the following abbreviated form: 


n 
NZ = XG + Ej Xj — Max 
j=l 
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Constraints: 
n 
>, Naijx; <b,+6,; ; t=1,2,..,m 
j=1 
xj 20 
It is given in the following detailed form: 
Find 
NZ = Neyx1, + NepX2 +++: + Ney xX, — Max 
Constraints: 
NQ44X1 + NQ42X2 ++» + NQyyn Xn S ND, 
NQz4X1 + NAz2X2 + +++ NAgnXy S Nb 


NQmiX1 + NAm2xX2 +1 + NAmynXn S ND 
X4)X2).4,Xy 2 O 


Linear models can also be expressed using matrices, and 
therefore the neutrosophic linear model given in the canonical 
form can be written using matrices as follows: 


Find 
NZ = NC X — Max 
Constraints: 
NAXSNB 
X2>0 


re 
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Nay, Naj42...Nazn Nb, Nex 
NA= Nay, Naz ... NA2zn NB = Nb» NC = oe 
Nami Nan? --»NQmn Nby, 


2-1-3- The standard neutrosophic formula for the linear 
model: 


The standard form plays an important role in finding a solution 
to linear programming problems, as the issue of searching for a 
solution to a linear programming problem has _ been 
transformed into the process of searching for a solution to a set 
of linear equations consisting of n equations with n+m 
unknowns, and solving this sentence is useful if it is possible, 
i.e. If it fulfills the conditions of non-negativity x; 2 0 , then 
the optimal solution for the linear model is the ideal values of 
the variables that fulfill the constraints and give the objective 
function the greatest or smallest possible value according to the 
text of the problem being solved. The standard neutrosophic 
formula is given in the following abbreviated form: 


Find 
NZ= XG + Ej )x; — (Max or Min) 


Constraints: 


>, Nai; = bi t 5 ; b= 1,2); M 
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It is given in the following detailed form: 
Find 


NZ = NeyX1 + Ne2X2 + +++ + NCnX, — (Max or Min) 


Constraints: 


NQy4X4 + NQ42X2 fies NaQinXn => Nb, 
Naz4X4 + NQp2X2 five NQgnXn = Nb, 


NQmiX1 + NAm2X2 +1 + NAmnXn = ND 
Kip Ks weg hye O 


Linear models can also be expressed using matrices, and 
therefore the neutrosophic linear model given in the general 
form can be written using matrices as follows: 


Find: 
NZ = NC X — (Max or Min) 


Constraints: 


NAX=NB 
X20 
Where: 
Nay, Nay2...Nain Nb, Nc, x4 
NA= Nay, NQp2 ...NA2n NB = Nb, NC = Nc X= x2 
Nang Namnz NQinn Nb New — Len 


Here we note that all constraints are of the equality type, except 
for the non-negativity constraints, which remain inequalities. 
Also, the right side of each equality constraint must be non- 
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negative, and all decision variables must be non-negative. The 
objective function in the standard Neutrosophic form can be a 
maximization function or a minimization function. 


2-1-4- The symmetrical formula of the neutrosophic linear 
model: 


We say of a linear program that it is in the symmetrical form if 
all variables are constrained to be non-negative and if all 
constraints are given in the form of inequalities , the 
inequalities of the constraints of the maximization problem 
must be in the form (<) less than or equal to, while the 
inequalities of the constraints in the minimization problem 
must be In the form (=) is greater than or equal to, then we 
write the neutrosophic symmetric formula in one of the 
following two forms: 

First figure: 

The neutrosophic symmetric formula for the linear 
mathematical model is given in the abbreviated form as 
follows: 


n 
NZ = XG + &;)x; — Max 
j=l 


Constraints: 


It is given in the following detailed form: 
Find 
NZ = Ne,x1, + Ne2X2 +++: + Ney X, — Max 
230i 
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Constraints: 
NQy4X4 + NQ42X2 + ee + NainXn < Nb, 


NQp4X4 + NQp2X2 + he + NQpznXn = Nb, 


NQmiX1 + NAm2X2 ++ + NAmnXn = NDm 
Ny Noy ti, Xq 0 
Using matrices as follows: 
Find 
NZ = NC KX — Max 


Constraints: 


NAX<NB 
X20 
Where: 
Nay, Nay2...Nain Nb, Ncy x4 
ae Nay, Nag ...NQA2n NB = Nb» Nc Nc X= x2 
Nama Naz NQinn Nb Nn us 


Second form: 
The summary is as follows: 


The neutrosophic symmetric formula for the linear 
mathematical model is given in the abbreviated form as 
follows: 


n 
NZ = YG + Ej )x; — Min 
j=l 


Se 
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Constraints: 


n 
>, Naix; >b,+6; ; t=1,2,..,m 
j=1 
xj 20 
It is given in the following detailed form: 
Find 
NZ = NeyXxX1 + Ne2X2 +++: + NC,X, — Min 

Constraints: 

NQ44X1 + NQ42Xq ++» + NQy Xn = Nb, 

NA z4X1 + NAz2X2 + +++ + NAgnXy = Nbz 


NQmiX1 + NAm2xX2 +1 + NAmnXn = Nb 
Xj X95) cay ty, &0 
Using matrices as follows: 
Find 
NZ = NC X — Min 


Constraints: 


NAX=WNB 
X20 
Where: 
Nay, Nay2..Nayn Nb, Ncy x4 
NA= Nay, NQy2 ... NA2n NB = Nb, NC = Nc X= x2 
Nam Naqz NO Nbn Nen i 


a) 
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2-2- How to move from one formula to another: 
Some presentation of the formulas of the neutrosophic linear 
models. It should be noted that we can move from one formula 
to another by following the following elementary 
transformations: 

e Converting the minimum value of the objective 
function f(x) to a maximum value by multiplying it 
by (-1) we get (- (f(x)). 

e If the inequalities were of the form (greater than or 
equal to) they will be converted to the form (less than 
or equal to) by multiplying both sides by (-1), and 
vice versa. 

e The equality constraint can be converted into two 
inequalities of different direction. 

e If the left side of an (inequality) constraint is given in 
absolute value, it can be converted into two regular 
inequalities. 

e Constraint inequalities of the type (greater than or 
equal to) are converted to an equality constraint by 
subtracting an appropriate positive variable (e., 
artificial variable) from the left side of the inequality 
and this variable is entered into the objective function 
with zero coefficient. 

e Constraint inequalities of the type (less than or equal 
to) are converted into an equality constraint by 
adding an appropriate positive variable (i.e., slack 
variable) to the left-hand side of the inequality and 
then this variable is entered into the objective 
function with zero coefficient. 
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e If one of the decision variables x is not constrained 
by the non-negative condition (that is, it can be 
negative, positive or zero), then it can be expressed as 
the difference between two non-negative variables 
x',x as follows x = x’ —x'and x'’,x >0 

2-3- Examples of the above: 

The linear models in all examples are given in detailed form: 
Example 1: 

Call the following Neutrosophic Linear Programming in its 

general form: 


Min NL = (3 4 €,)x, — (3 £ €3)X_ + (7 £ €3) x3 
Constraints: 


Xy +X. + 3x, 54046, 
X1 + 9x2 — 7x3 = 504 62 
5X, + 3x, = 204+ 63 
[5x2 + 8x3| < 100+ 6, 
X1,X, 20 
where €; It is indeterminate and could be 
€&j € [Ayj,Azj| or & € (Ay; Aaj} sf = 1,2,3. 
Also, the values that express the available possibilities 6; are 
neutrosophic values. This means that 
It is indeterminate and could be 
6; © [Mii Mai] or 6; € {Hi Mai} st = 1,2,3,4 
To convert the above problem into the neutrosophic canonical 
form, we perform the following transformations: 
¢ The objective function is a function of minimization that we 
turn into a function of maximization 


Min NL = (34+ €,)x, —(3 4 €)x. + (7 £ €3) x3 


“OAs 
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transformed into 
Max NZ = —(3 + €:)x, + (3 £ €))x2 — (7  €3) x3 
¢ The second constraint is given (greater than or equal to) is 
converted into (less than or equal) by multiplying both sides by 
(-1) we get —x, — 9x2 + 7x3 < —(50 + 62) 
¢ Third constraint 5x, + 3x, = 20 +63 transformed into two 
entries 
5X1 + 3x2 < 204 63 
5X1 + 3x2 = 204 63 
Then we turn the constraint 5x, +3x, =20+63 into 
—5x, — 3x2 < —(20 + 63) 
* The constraint [5x2 + 8x3| < 100 + 6, is equivalent to the 
two inequalities 
5x2 + 8x3 < 1004+ 6, 
—5xX2 — 8x3 < 100+ 6, 
¢ The variable x3 is not restricted by the non-negative 
constraint, so it is replaced by the following assumption 
X3 = X4 — X3wherex4, x3 > 0. 
The canonical neutrosophic form becomes: 
Max NZ =—G+¢e,)x,+ GB te)x. —(7 + €3)(X3 — x3) 
Constraints: 
x, + x2 + 3(x} — x3) $4046, 
—x, — 9x5 + 7(x5 — x3) < —(50 + 6,) 
5x1 + 3x2 < 20463 
—5x, — 3x2 < —(20 + 63) 
5x, + 8(x3 —x3) < 100+ 6, 
—5x, — 8(x3 —x3) < 100+ 6, 


I W 
X1,X2,X3,X3 2 0 
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Example 2: 


A factory produces four types of products Sj, Sz,,S3,S4. For 
this purpose, the following raw materials are used M,, M2, M3 . 
The factory management wants to study the optimal 
organization of production during a period of time (for 
example, a month) and determine the monthly production for 
each product in order to achieve a maximum profit, bearing in 
mind that the profit is directly proportional to the number of 
units sold of the products. The available quantities of raw 
materials needed for each product and the profit have been 
showed in the following table: 


Products Product Type 
Matesals Available 
Sy S2 S3 Sy Quantities 
M, 1.5 1 2.4 1 3000+ 6, 
Mz 1 5 1 3.5 9000 + 6, 
M3 1.5 3 3.5 1 7000 + 63 
win one 
4+é€, | Bt& | 5t€ |) 648 
product 


Let us suppose x, is the number of units produced from the 
first type, x2 is the number of units produced from the second 
type, x3 is the number of units produced from the third type, 
and x, is the number of units produced from the fourth type 
during the production period (a month for example), and 
accordingly, the consumed quantity of the raw material M, in 
the production of the four varieties will be: 


1.5x, + X2 + 2.4x%3 + x4 
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And it must not exceed 3000+ 6, from the available quantity, 
that is: 


1.5%, +x + 2.4%, +x,< 300046, (1) 


Likewise, the amount of raw material M> consumed in the 
production of the four types is: 


X1 + 5X2 +X3+3.5x, < 9000 + 6, (2) 


And the amount consumed of the raw material M; in the 
production of the four types is: 


1.5x, + 3x2 + 3.5x3 + x4 < 7000 + 63 (3) 
In addition, the produced quantities must be non-negative, 1.e.: 
X4,%X2,X3,X, 20 (A) 
And they are what is called non-negative conditions. 


Thus, we have identified all the constraints imposed on the 
variables of the problem. 

We now define the objective function. If quantified units 
X4,X2,X3,X4 Of species are produced in order, then the profit 
during the productive period will be: 


NZ = (4 + E1)X4 + (8 ke Ey) Xo + (5 + E3)X3 + (6 BE E4)X4 
It represents the target function. Therefore, the mathematical 


model of the problem is: 
Max NZ = (4+ €)x, + (8 £ €&)x. + (5 + €3)xX3 + (6 + &4)Xy4 


Constraints: 
1.5%, +x, + 2.4x, + x, < 3000+ 6, 
X1 + 5x2 + x3 +3.5x, < 9000 + 6, 
1.5x, + 3x2 +3.5x3 +x, < 7000 + 63 
X4,X2,X3,X, =O 
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We have obtained a neutrosophical canonical linear model 
using the appropriate transformations, which can be written in 
the following neutrosophical standard form: 


Max NZ = (44 €,)x, + (8+ €&)x2 + (5 + €3) x3 
+ (6 + €3)x, + Oy, + Oy, + Oy; 


Constraints: 


1.5%, +%X%2+2.4x%3 +x, + y, = 300046, 
X1, + 5x2 +X3 + 3.5x, + yp = 9000 + 6, 
1.5xX; + 3X5 + 3.5X%3 + X4+ yo = 7000 + 6, 
X14, X2,X3,X4, =O 


Conclusion: 


The indeterminacy that we added to the data described by the 
linear model provides us with neutrosophical linear models that 
simulate reality and take into account most of the changes that 
could occur in the operating environment of the system 
represented by the linear mathematical model, and enable us to 
continue studying linear programming topics such as finding 
accompanying programs that need to be developed. The 
mathematical model in the symmetrical form, solving linear 
models using the simplex method that requires developing 
models in the standard form, and other linear programming 
topics. 
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Chapter ITI 
The graphical method for finding the optimal 
solution for neutrosophic linear models 


Introduction 
3.1. Graphical method for solving linear models 


3.2. Graphical method for finding the optimal solution for 
neutrosophic linear models 


3.3. Non-negative constraints for optimal solution of some 
neutrosophic linear models using the graphical method 


3.4. Neutrosophic linear mathematical model conclusion 
Conclusion. 
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Chapter ITI 
The graphical method for finding the optimal 
solution for neutrosophic linear models 


Introduction 

After introducing the linear models and their different formulas 
based on the concepts of neutrosophic science, in this chapter 
we present the neutrosophic graphical method that we use to 
solve the neutrosophic linear models. 

The graphical method represents the model graphically and is 
one of the simplest methods for solving linear programming 
problems. However, it is not sufficient to solve all linear 
programming problems, as linear programming problems often 
contain a large number of variables, and the use of the 
graphical method is limited to the following cases: 

e The number of unknowns is n= 1, orn=2, orn=3. 

e In linear models whose constraints are equal constraints, if 
the number of unknowns and the number of equations meet 
one of the following conditions: n—m= 1 orn—m=2or 
n-m=3. 

Here we can transform the model into a function of one 
variable, or two variables, or three variables; by using the non- 
negative constraints that the variables of the linear model have. 
In this research, we present a reformulation of the graphical 
method for solving linear models using neutrosophics, as well 
as the graphical method for solving linear models where the 
constraints are equal, and the difference between the number of 
unknowns and the number of constraints is equal to one, two, 
or three. 
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3.1. Graphical method for solving linear models 
We find the optimal solution by following the steps below: 


1. We determine the half-planes defined by the inequalities 
of the constraints by drawing the straight lines resulting 
from the transformation of the inequalities of the 
constraints. To do this, we specify two points that fulfill 
the constraint, and connect the two points to obtain the 
straight line that corresponds to the constraint. This 
straight line divides the plane into two halves to 
determine the half-plane that satisfies the constraint. We 
select a point at the top of the mapping from one of the 
two half-planes. We substitute the coordinates of this 
point into the inequality. If it is satisfied, then the region 
in which this point is located is the solution region. If it 
is not satisfied, then the opposite region is the solution 
region. 

2. We define the common solution region, i.e., the region 
resulting from the intersection of the halves of the 
planes defined by constraint inequalities. This region 
must be non-empty so that we can proceed with the 
solution. 

3. To represent the objective function, we note that its 
relation contains three unknowns, Z, X41 Xp. 
Therefore, we need to know a value for Z that is 
unknown to us. Here we assume a value, let it be Z, = 
0, draw the equation of the objective function Z, specify 
another value, let it be Zz, and represent the equation. 
We get a line that is parallel to the first line, and if we 
continue like this, we get a series of parallel lines that 
represent the target function. 
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4, 


3.2; 


a: ~ “EE 
We draw ray C = | | where c, is coefficient of x, and 


Cz is coefficient of x2 in the objective function 
statement, and the direction of its increasing function is 


4 C 
the direction of ray C = | Bell and the direction of its 


decreasing function is the opposite direction. This ray, 
ie., the drawing is done according to the type of 
objective function (maximization or minimization). To 
put it more clearly, we find the optimal solution point by 
drawing the line representing Z, parallel to itself 


> yc 
towards the ray C = | ie to find the maximum value of 


the objective function (and reversing this direction to 
find the smallest value), until it passes through the last 
point of the common solution region and this point is the 
optimal solution point, which is located at the 
boundaries of the common solution region and any other 
displacement, no matter how small, takes it out of it. 


Graphical method for finding the optimal 
solution for neutrosophic linear models 


From the definition of neutrosophic linear models, we find that 
we can apply the same previous steps to obtain the optimal 
solution, which is a neutrosophic value suitable for all 


conditions. We illustrate the above with the following example: 


Example 1 


A company produces two types of products A, Az and uses 


three 


types of raw materials B, B,,B3in the production 


process; the available quantities of each of the raw materials 


are B; ;t = 1,2,3, the quantity required to produce one unit of 
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each products A; ; j = 1,2, and the profit derived from one 
unit of each of the products A, Az is shown in the following 
table: 


products A, rm available 

raw materials quantities 

B, 6 4 36 

B2 Z 3 12 

B3 5 0 10 
profit [6, 8] [2,4] 

Table Issue data 

Requirement 


Determine the quantities that must be produced of each product 


A; ; Jj =1,2 in order for the company to achieve maximum 


profit: 
Solution 


Suppose x; is the quantity produced from the product, 
where j = 1,2, then we can formulate the following 
neutrosophic linear mathematical model: 


Z = [6,8]x, + [2,4|x. ~ Max 
Constraints 
6x, + 4x, < 36 (1) 
2X, +3x2<12 (2) 
5x, <15 
X1,X, 20 
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The previous model is a linear neutrosophic model because the 
coefficients of the variables in the objective function are 
undetermined. To find the optimal solution for the previous 
model, we will apply the graphical method according to the 
following steps: 


The first constraint 
We draw the straight line representing the first constraint: 
6x, + 4x2 = 36 
We impose: 
x, =0>5> 4x, =365x%,=9 
We get the first point: A(0,9). 
We impose: 
Xp =0 > 6x, = 36>x, =6 
We get the second point: B(6,0) 


We take a point at the top of the designation from one of the 
two halves of the resulting plane after having drawn the 
straight line through the two points A(0,9) and B(6,0). Let it 
be the point O(0,0) and substitute it in the inequality of the 
first entry. We find that the inequality is satisfied 1.e., the half 
of the plane to which the point O(0,0) belongs is half of the 
solution plane of the first-constraint inequality. 


We proceed in the same way for the second and _ third 
constraints and obtain the following graphical representation: 
Figure No. (1) 
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Figure No. (1) Graphic representation of the limitations of the linear 
model in Example 1 


After we have shown the constraints, we notice that the 
common solution area is bounded by the polygon whose 
vertices are the points, (0,0) , E(3,0) , M and C(0,4) 


The point M is the intersection point, and the second and third 
constraints coordinates are obtained by solving the following 
two equations: 


2X4 + 3X2 = 12 
5x, = 15 
We find : M(3,2) 


Substituting the coordinates of the vertex points into the 
objective function expression, we get: 


Zo =0 
Ze € [12,16] 
Zu € [22,32] 
Zc € [8,16] 
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This means that the highest value of function Z is reached at 
point (3,2),i.e., the company must produce three units of the 
first product and two units of the second product, then it will 
achieve the maximum profit. 


Max Z =Zy € [22,32] 
Note 


The process of substitution of the objective function by the 
coordinates of the points of the vertices of the common 
solution area is possible when the number of points is small, 
since we can easily substitute them in the objective function, 
and the point that gives the best value for the objective function 
is the optimal solution, but when there is a large number of 
constraints, we get a large number from the vertical points 
located on the perimeter of the common solution region. In this 
case, the method of determining the coordinates of all these 
points and substituting them into the objective function 
becomes impractical. Therefore, we resort to the representation 
of the objective function and the determination of the optimal 
solution point as we have already mentioned. 


3.3. Non-negative constraints for optimal solution of 
some neutrosophic linear models using the 
graphical method 


Example 2 


Find the optimal solution for the following linear neutrosophic 
model: 


Z =X, — Xp — 3X3 +X, + ([2,5]x5 — x6 + 2x, — [10,15] — Max 
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Constraints 
Xy—Xz+xX3=5 (1) 
2X4 —X_2—-X3—-X,=-11 (2) 
Xy +X2—-—X, = —4 (3) 
Xz +x, =6 (4) 


2X4 = 3X2 = X6 + 2X7 = 8 (5) 
X41, X2,X3,X4,X5,XE,X7 =O 
Solution 


We note that the number of constraints is m=5 and the 
number of variables is n = 7, which means that n — m = 2. 
Therefore, based—on the non-negative constraints, find the 
optimal solution for the previous model using the graphical 
method according to the following steps: 


1- We calculate five variables in terms of only two 
variables. 

2- Since the variables of the linear model satisfy the non- 
negative constraints, then we obtain from the variables 
that we calculated five inequalities of the type greater 
than or equal to. 

3- The objective function with only two variables is 
obtained by substituting the five variables. 

4- We write the new model, which is a linear model with 
two variables, so that the optimal solution can be found 
graphically. 


We apply the previous steps to Example 2. 
We find: 
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X, =5—-x, +x, (1)’ 
Xq = 3x, -2x,+6 (2)! 
Xe =X, +x, +4 (3)' 
Xe = 6— Xz (4)' 


Xy = 7— xX, +X, (5)' 
Substituting in the objective function, we get: 
Z = [1,4]x, + [3,6]x2 + [8,25] 
Since X3,%4,%s5,X6,X7 =O from (1)’, (2)’, (3)’, (4)’, (5)’, we 
get the following set of constraints: 
5—X, +X, 20 
= 32) lx 3 2 0 
Xy+%X%2+420 
6—Xx, 20 
7 — Xy +X2 20 
Neutrosophic linear mathematical model 
Requirement 
Z = [1,4|x, + [3,6]x, + [8,25] > Max 
Constraints 
5 —X, +X, 20 
3X, — 2x, +620 
Xy+xX2+420 
6—Xx, 20 


7 — X% +X S| 0 
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X1,X2 20 


The model has two variables, so the optimal solution can be 
found graphically by following the same steps indicated in 
Example (1). 


The required graphic representation is found in Figure No. (2). 


A we 


Figure No. (2): Graphical representation of the constraints of the 
linear model in Example 2 


Region D is the region of joint solutions and is defined by the 
polygon OBRSC, where,O(0,0) , (5,0) , C(0,3) , and for the 
two points R ,S we find: 


Point R is the point of intersection of the first and fourth 
entries. 


We obtain its coordinates by solving the set of equations: 
5— xX, +x, =0 
6—x,=0 
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We get: R(11,6) 


Point S is the point of intersection of the second and fourth 
entries. 


We obtain its coordinates by solving the set of equations: 
3X, — 2X, +6=0 
6—x,=0 
We get: S(2,6) 


Since the optimal solution is located at one of the vertices of 
the common solution region, we substitute the coordinates of 
these points with the objective function: 


At point 0(0,0) 
Zo =0 

At point B(5,0) 

Zp = [13,45] 
At poin tR(11,6) 

Zp € [37,105] 
At point $(2,6) 

Z; € [28,69] 
At point C (0,3) 

Zc € [17,43] 


The greatest value of the objective function is at the point 
R(11,6) that is x, = 11 andx, = 6. 


We calculate the values of the remaining variables by (1)’, (2)’, 


(3), (4), GY. 
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We find: x3 = 0 ,x, = 27 ,x5 = 21,%, =0,x, =2. 


Substituting in the objective function of the original model we 


obtain the maximum value of the Z function, which is: 


MaxZ € [68,126] 


Important Notes 


|- 


A vertical point in space R” is covered by the graphical 
solution. The ideal solution pertains to a vertical point, 
which is the outcome of several lines or planes 
intersecting, therefore the number of non-existent 
components is at least nm — m components. 

Certain conditions that are irrelevant to the solution 
process might be included in the model. 

When one of the sides of the common solution area that 
passes through the ideal solution point is parallel to the 
straight-line Z=0, the ideal solution can be a single point 
or an infinite number of points. Thus, when the objective 
function is represented by a straight line, this line will 
apply to the parallel side, and all of the infinitely many 
points on that side will be perfect solutions. 

We say that the objective function has an endless 
number of acceptable solutions that offer us greater 
values of Z if the region of acceptable solutions is open 
in terms of growing the function Z, meaning that we 
cannot stop at a particular perfect solution. 

The situation in which there is no perfect (acceptable) 
solution and the zone of alternatives is an empty set (the 
problem is impossible to solve) when the conditions 
conflict. 
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Conclusion 


This chapter covered both the graphical method and a method 
that is rarely covered in references on classical operations 
research: using non-negative constraints to graphically find the 
optimal solution for some neutrosophic linear models. 
However, it should be noted that there are some neutrosophic 
linear models that have two variables. In these cases, it may be 
difficult to reach the common solution region or to determine 
the optimal solution once the common solution has been found, 
so it is preferable to use the simplex neutrosophic method. The 
researcher must choose the best approach for the model he 
wishes to solve because the main objective is to arrive at the 
best solution. 
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Chapter IV 


The simplex direct neutrosophic algorithm for 
finding the optimal solution for linear models 


Introduction. 


4-|- The neutrosophic linear models set in the symmetrical 
form and of the Max type. 


4-2- The neutrosophic linear models are in symmetric form 
and are of type Min. 


Conclusion. 
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Chapter IV 


The simplex direct neutrosophic algorithm for 
finding the optimal solution for linear models 


Introduction: 


Linear programming is the method that helps in selecting 
decisions and approving the best program for independent 
activities, taking into account the available resources. Linear 
programming is used in solving problems in which the goal is 
specific, such as securing a maximum profit, securing a 
minimum cost, or saving the greatest time or effort ... Etc., 
noting that the problem of linear programming consisting of a 
linear function and knowledge of a set of inequalities or 
equations (constraints) is characterized by the presence of a 
large number of acceptable non-negative solutions, and what is 
required is to find the optimal solution from among these 
solutions. To reach this solution, we rely on the information 
that we presented it in the first chapter when we studied non- 
negative solutions to the system of neutrosophic linear 
equations. Then we used the simplex method, which is the 
mathematical basis for the direct simplex algorithm used to 
find the optimal solution for the linear models that we will 
present in this chapter: 


Direct simplex algorithm for solving neutrosophic linear 
models: 
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The direct simplex algorithm consists of three stages: 


a- The stage of converting the imposed model into an 
equivalent systematic form. 

b- The stage of converting the regular form into a basic 
form to obtain the non-negative basic solutions. 

c- The stage of searching for the optimal solution required 
from among the non-negative basic solutions. 


In this chapter, we will use the direct simplex algorithm to 
find the optimal solution for the neutrosophic linear models, 
which were presented in the second chapter of this book, 
and here we distinguish the following cases: 


1- The neutrosophic linear models are in symmetric form 
and are of type Max. 
2- The neutrosophic linear models are in symmetric form 
and are of type Min. 
3- The neutrosophic linear models are given in the general 
form. 
Using the direct simplex method to find the optimal 
solution 


4-1- The neutrosophic linear models set in the 
symmetrical form and of the Max type: 


The neutrosophic linear model of the Max type is written in the 
symmetrical form, as we mentioned in the second chapter, in 
the following form: 


Find: 
NZ = Neyx1, + Ne2X2 +++: + NCyX, — Max 
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Within restrictions: 
NQy4X4 + NQ42X2 + oe + NainXn < Nb, 


Nap4X4 + NQp2X2 + is + NQznXn < Nb, 


NQmiX1 + NAm2xX2 t+ NAmynXn S ND 
Ky XopanyXy 0 
The study is carried out according to the following steps: 


1- We write the model in standard form; we get the 
following form: 


Find 
NZ = NeyX1, + NepX2q +°1° + NCyXy + 0.1 + 0. V2 +++ +0. Yn — Max 
Within restrictions: 
Nay1X1, + Nay2X2 + °° + NainX, + ¥, = Nb, 
NQz14X1 + NQg2X2 +9 + NAgnXn + V2 = Nbz 


N@piX{ TNO yoxXs 1: + NG Xp + Vm = ND 


Ke Royce Ky Vin Vir Vn © 


2- We convert the model to the basic form. Here we notice 
that the additional variables serve as an initial base from 
which we start searching for the optimal solution. 
Therefore, we arrange the model information in the 
following table: 
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oe x4 X2 x Available 
basicw. | 41 | 42° | «ee - Vy | V1 | | Vm ame 
V4 Nay, Nay2 ara Nain 1 0 arate @) Nb, 
V2 Na, Nay Tr NQzn 0 1 aa @) Nb, 
Yn Nami | NQAm2| -- |NQmn| 0} O0]..) 7 ND 
objective 
function Ne, | Nez . | Ne, |O}0;0] 0 | Z-O 


Table No. (1): Basic information of the model 


We have a first base consisting of the variables y,, y2,.--,Vm; 
then the variables x1, X2,...,X,are non-base variables and we 
move to the next step: 


3- We determine the appropriate variable from the 
equations and insert it into the rule by studying examples 
of the variables in the row of the objective function Z. 
Since the objective function is a maximization function, 
we choose the largest positive values in the row of the 
objective function. In other words, we take: 


Max(Nc,, NCo, ...,NCy) = NC, 
For example, let it be Nc, corresponding to the variable x,. 


Thus, we have determined the pivot column. This means that 
the variable x, will enter the base to determine the variable that 
will exit from the base, and therefore the pivot line. We 
calculate the following indicator: 

= > 0; Na;; > 0,Nb; > 0 
Nais Nats 


6 € Min| 
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The element located at the intersection of the pivot column 
with the pivot row is the pivot element. 
- We divide the pivot row by the pivot element, we get: 
Na, Nay Nays_4 Naes41 Nagy Nb, 
Na,, Na; Na ’ ’ Nas’ Nas’ Nats 
- We make all the elements of the pivot column equal 
to zeros, except for the pivot element, which is equal 


to one. 
- We perform the appropriate calculations to calculate 
the current of the new table using the following 


relationships: 
Najs _ NajjNa¢s = Na,jNdis 


Na;; = Na;; — Na,; —— = 
J tj 
? Nats Nats 
; Na,;, Nb,Na,, — Nb,Na; 
Nb, = Nb; — Nb, == 
Nats Nas 
; Na  NcjNas — Nc,Naz; 
Nc; = Nc; — Ne, = 
Nats Nats 
We get the following table: 
riables : 
es s a ; Available 
; . ee - _ nm | ¥1 | ¥1 | | Ym quantities 
basic 
Vi Nai, | Nays | «| NQyo4 | Nain | 1 | O |. | 0 Nb 
V2 Naz, Nay2 ane Nabs-1 wee Nazn 0 1 ane 0 Nb; 


Ym Nam | NAmz2 | «. NGipsa4 «| NQmn | 0 0 |... 1 Nb}, 
objective , ; ; 
function Ney | Neg | | Nes-1 «| Ne, | 0 | 0 | 0} O NZ 


Table No. (2) The first step in the simplex direct neutrosophic 
algorithm 
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- We apply the stopping criterion of the Simplex algorithm 
to the objective function row in Table No. (2). 


Stopping criterion: 


Since the objective function is of the maximize type, the 
objective function row in the table must not contain any 
positive value, (but if the objective function is of the 
minimization type, the objective function row in the new table 
must not contain any negative value), in If the criterion is not 
met, we return to step No. (3) and repeat the same steps until 
the stopping criterion is met and we obtain the required optimal 
solution. Thus, we obtain new non-negative neutrosophic basic 
solutions and non-basic (free) solutions equal to zero. The ideal 
solution is written in the following form: 


(Nb, Nb3, ..., Nb, 0,0, ...,0) 


The following table represents the final solution if the basic 


solutions are: (x1, Xz, «.,Xm) 
variables Available 
xX, |x x x 
Ea its” i le an sa ¥2 Ym quantities 
x4 1] 0 0 | Naimss Nain | NB | NBit NBim Nbi 
Xp 0/1 0 | Nagm+1 Nagn | NB21 | NB22 NBom Nb} 
0 0 
Xm 0 0 1 N@nm+1 NQmn NBm1 NB m2 NBimm Nb», 
objective | 0 | 0 Oo). NG ae Ney, | Nay | Nao Nam NZ' 
function 


Table No. (3) The final solution in the simplex direct neutrosophic 
algorithm 


Where Nf;; and Nq, are the examples of the additional 


variables in the constraints and in the objective function after 
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performing the aforementioned iterative operations, the optimal 
solution is 


Xx, = Nbj, xX. = Nbj, ..,Xm = Nbm 
Which gives the maximum value of the following objective 
function: 
NZ' = Nc, Nb, + Nc2Nb, +++ + NcyNbi, 
We explain the above through the following example: 
Examplel: 
Problem Statements Classical Values: 


A company produces two types of products A, B using four raw 
materials F,, Fj, F3,F,. The quantities needed from each of 
these materials to produce one unit of each of the two 
producers A,B, the available quantities of the raw materials, 
and the profit returned from one unit of both products are 
shown in the following table: 


Products | Required quantity per Available 
unit quantities of 

Raw Materials . x a 
F, 2 3 19 
Fy 2 1 13 
F3 0 3 15 
F, 3 0 18 

Profit Returned per unit 7 5 


Table No. (4) Classic data for the issue 
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Required: 


Finding the optimal production plan that makes the company's 
profit from the producers A, B as large as possible. 


We symbolize the quantities produced from the product A with 
the symbol x,, and from the product B with the symbol x2 , 
after building the appropriate mathematical model and solving 
it, we conclude that x, = 5,x, = 3 , and hence the maximum 
profit MaxZ = 50 of monetary unit. 


Problem Statements neutrosophic Values: 

A company produces two types of products A, B using four raw 
materials F,, F,,F3,F,. The quantities needed from each of 
these materials to produce one unit of each of the two 
producers A,B, the available quantities of the raw materials, 
and the profit returned from one unit of both products are 
shown in the following table: 


Products | Required quantity per unit Available 
A B quantities of the 
raw materials 
Raw Materials 
un 2 3 [14,20] 
Fo 2 if [10,16] 
Fs 0 3 [12,18] 
F4 3 0 [15,21] 
Profit Returned per [5,8] [3,6] 

unit 


Table No. (5) Neutrosophic data for the issue 
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Required: 


Finding the optimal production plan that makes the company's 
profit from the producers A, B as large as possible. 
Symbolize the quantities produced from the product A with the 
symbol x,, and from the product B with the symbol x2 , the 
problem will be reformulated from the neutrosophic 
perspective as follow: 
NZ = [5,8] x, + [3,6] x. — Max 
Within restrictions: 
2x, + 3x2 + y, = [14,20] 
2x, +x2+y, = [10,16] 
3x2 + y3 = [12,18] 
3x, + v4 = [15,21] 
X1,X, = 0 
The above program needs to reformulated to an equivalent 
form by adding slack variables: 
NZ = [5,8] x, + [3,6] x. + Oy, + Oy, + Oy; + Oy, — Max 
Within restrictions: 
2x, + 3x2. + y, = [14,20] 
2x, + x2 + y2 = [10,16] 
3x2 + y3 = [12,18] 
3x, + y¥, = [15,21] 


X41%X2) V1» V2» V3 V4 = 0 


We arrange the previous information in the following table: 
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x4 Vi Ve | Voe| Va Available 

quantities 
Vi 3 it 0 0 0 [14,20] 
V, 1 0 1 0 0 [10,16] 
[12,18] 


objective 
function 


Table No. (6): The first step in the simplex method 


We note that the additional variables form an initial base 
consisting of the variables (y,,y,,y3,y,). Then we consider the 
variables (x,,x2) are non-basic variables and we move to the 
next step: 


We determine the appropriate variable from the equations and 
insert it into the rule by studying the examples of the variables 
included in the expression for the objective function NZ. Since 
the objective function is a maximization function, we choose 
the variable with the largest positive examples from the last 
row in the table, that is, from the row of the objective function. 
In other words, we take 


Max([5,8], [3,6]) = [5,8] 


It is clear that versus to the column of x; , meaning that the 
variable x, should be placed instead of one of the basic 
variables. Now to demonstrate which basic variables should be 
ejected, the following calculation has been done: 

[14,20] [10,16] [15,21] [15,21] 


Oe ae 
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The value of @ indicates that the row versus to the variable y4, 
and the element positioned in the cross row/column is 3 which 
is the pivot element, divide the elements of the row versus to 
y4 yields: 

30000 1 [15,21] 

=) sisi sis is) —  LD || 

3a 3 3 3 9 3 
Then we make all the elements of the pivot column equal to 
zero, except for the pivot element, which is equal to one. We 
perform the appropriate calculations using the following 
relationships: 


Najs NajjNa¢s = Na,jNadis 


ro _ 
peti nat = a ae = Nats 
; Na, Nb,Na,. — Nb,Na, 
Nb, = Nb; — Nb, 2 = 
Nas Nas 
Na,; Nc,;Na,. — Nc.Na;; 
Nc = Nc; — Nc, a Paraiso 3 
Nas Nats 


We get the following table: 


ariables Available 
Gace ay Yi | Yo | Ys V4 quantities 


=Z 

vi 0 1 0 0 a [4,6] 
=v 

Y> 0 0 1 0 = [0,4] 


1 1 0 | 0 | 0 = [5,7] 
objective —8 -5 | NZ 
0 0 0 ) ——— 
function 3’ 3 1) - [25,56] 


For No. (7), the second step is the simplex method 
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It still there is a non-negative value in the row of the objective 
function (i.e., [3,6]) which is versus to the xz column, this 
leads to the fact that the variable xz should be entered into the 
basic variables. Now, the question is: which basic variable 
should be ejected? Track the following calculation to answer 
this question, 

[4,6] [0,4] [12,18] [4,6] 4 


6 € Min |—— ,——_ = | 
= ie oe ee 3 134! 


which is versus to the slack variable y,, the pivot element equal 
3, hence the row versus to y, should be divided by 3 , the 
required calculations yield the following tableau: 


Variables Available 
haste a4 x2 V4 Y2 Y3 V4 quantities 
1 —2 
X2 0 1 — 0 0 —- fe 
3 9 ) 
i —4 - 
0 0 = 1 0 — — 
2 
Vs 0 | 0 -1 0 1 5 [8,12] 
1 
a6 1 0 0 0 0 5 [5,7] 
objective —6 NZ 
0 0 -2,-1 0 0 —,- 
function a Wy - [29,68] 


Table No. (8) Final solution 


It is clear from the row of the objective function that all the 
elements are either zero or neutrosophic negative numbers, this 
means that we have reached to the optimal solution is: 


4 4 
xt € [5,7],x3 € |=. 2],93 C 2,93 € [8,12], yi = yi =0 


- 65 - 


Neutrosophic linear models and algorithms to find their optimal solution 


Substitute the above optimal solution into the objective 
maximum function, yields: 


MaxNZ € [5,8].{5,7] + [3,6]. B 2| = [25,56] + [4,12] = [29,68] 


which is identical to the result in the previous tableau. 


Substituting the optimal solution into the constraints we find: 


2[5,7] +3 -. | +0 = [14,20] 


ais71+ £2] + [22] = Ho261 


3|- 2| + [8,12] = [12,18] 
3[5,7] + 0 = [15,21] 


We observe that the optimal solution satisfies all constraints. 


We summarize the previous results in the following table: 


results 


issue data 


Cy C2 b, bz b3 by x4 x2 


19 13 


issue data results 


Cin | Con bin bon b3n Dan Xin | X2Nn 


[5,8 | [3,6 | [14,20 | [10,16 | [12,18 | [15,21 | [5,7 
] ] ] ] | ] ] 


Table No. (9) Comparison between the results of solving the problem, 
classical data, and neutrosophical data 
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4-2- The neutrosophic linear models are in 
symmetric form and are of type Min. 


The neutrosophic linear model of the Min type is written in the 
symmetrical form, as we mentioned in the second chapter, in 
the following form: 


Find 
NL = Nceyx1, + Ne2X2 +++» + NCypX, — Min 
Within restrictions: 
NQ44X1 + NQ42X2 ++» + NA Xn = ND, 
NQz4X1 + NAg2X2 + +++ NAgnXy = Nb 


NQy1%1 + NQm2X2 + Shr + NGynXn = Nb», 
Kj Kopp Xycee 0 


We can search for the optimal solution by following one of 
the following methods: 


- The optimal solution for the previous linear model 
can be found by converting the objective function to a 
function of the maximization type, as we mentioned 
in the second chapter, by multiplying the line of the 
objective function by (-1), then we write the model in 
the standard form. Here we notice that there is no 
ready-made initial rule because the variables All 
additional ones are preceded by a minus sign, 
meaning we must first search for an initial solution 
and then improve this solution until we reach the 
optimal solution by following the same previous 
steps. 
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- Wecan also find the accompanying model, and it will 
certainly be like a symmetry of the maximization 
type, and then we find the optimal solution for it as 
we did previously, or by using the dual algorithm to 
solve the model and the accompanying model, which 
we will present in the seventh chapter of this book. 

- In such models, it is preferable to use the synthetic 
simplex algorithm, which will be presented in 
Chapter 6 of this book. 

- Also, we can find the solution without making any 
change in the objective function, but we make one 
modification to the aforementioned steps, which is 
that when we want to determine the anchor column, 
we choose the most negative element, and its column 
is the anchor column, and we follow the solution as 
we mentioned above, and the stopping criterion here 
is that all The line elements of the target function are 
either positive or zero. 

1- Neutrosophic linear models are given in the general form: 


The neutrosophic linear model is written in the following 
general form: 


Find 
NZ = Ney xX, + NepX2 + +++ + NCpnX, — (Max or Min) 
Within restrictions 
> 


Nady X1 + NAj2X2 +++ + NAinXn (=) Nb, it =1,2,...,m 


X4,Xp,++,Xyn 2 O 
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In a first step, we write the model in the standard form, and to 
write it in the basic form, we rarely have additional variables. 
Here we notice that some of the additional variables are 
suitable to be a basic variable, and others are not. Also, if there 
are some constraints of the type of equality, then there are no 
corresponding additional variables. As a result, there are no 
basic variables. A ready-made prototype, and we need to build 
an initial base from which to start searching for the optimal 
solution. Also, here it is preferable to use a simplex with an 
artificial base that will be presented in the sixth chapter of this 
book. 


Important Notes: 


In the final table for the model of the maximization type, if 
some of the examples corresponding to the free variables in the 
target function line are positive, this means that we have not 
reached the required ideal solution, and in this case, we must 
delete the free variables corresponding to the positive value. 
We return to step (2) and do the necessary. We repeat this. 
Operations until we get a target function line that contains 
nothing but zeros or negative (positive) numbers. Or we reach 
one of the following cases: 


a- There is no ideal solution because the solution region is 
open in the direction of increasing NZ, and we infer this 
from the absence of a positive element in the fulcrum. 

b- There is an infinite number of optimal solutions because 
the levels of the objective function NZ are parallel to one 
of the sides or surfaces of the common solution region. 
We infer this from the appearance of a zero in the last 
row corresponding to one of the free variables in the 
table of the last optimal solution. Then we can obtain 
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another optimal solution by changing the variable. with 
one of the basic variables, we will get another basic 
solution as a result. 

If there is no optimal solution, this happens because the 
constraints conflict with each other. We infer this from 
the absence of any positive element except for the 
constant Nb, in one of the lines. This means that the right 
side takes a negative value while the left side takes a 
positive value, and this only happens. When restrictions 
conflict. 

After obtaining the optimal solution, we must verify that 
it satisfies all the imposed constraints and gives us the 
same value for the objective function, by substituting the 
objective function and the constraints. 
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Conclusion: 


We conclude from the previous study and the results shown in 
Table (9) that when using neutrosophic data we obtain areas of 
any indeterminate values, and this indeterminacy is more 
accurate, simulates reality, and takes into account most of the 
changes that may occur in the operating environment of the 
system represented by the linear mathematical model, while the 
values that We obtain when solving according to classical data, 
they are specific values and do not take into account the change 
that may occur in the operating environment of the system 
represented by the mathematical model. Therefore, 
neutrosophical data provide us with a more general and 
comprehensive study than the study using known classical 
data, i.e. working using classical data It is no longer sufficient 
at the present time, because the development of science and the 
instability in the status of the facility’s work environment has 
placed before us a large number of cases that need quick and 
accurate treatment to avoid losses that the facility may be 
exposed to, which cannot be treated using classical data, and 
here comes the role of data. Neutrosophy, which provides us 
with more comprehensiveness in interpreting the results and 
helps us in obtaining the required accurate results. On the one 
hand, we focus on the necessity of choosing the appropriate 
algorithm to solve the model under study from among the 
algorithms that we will present in this book to save effort and 
time in searching for the optimal solution. 
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Chapter V 
Modified Neutrosophic Simplex algorithm to 
find the optimal solution for linear models 
Introduction: 
5-1- Steps of the modified simplex neutrosophic algorithm: 


Conclusion and results: 
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Chapter V 


Modified Neutrosophic Simplex algorithm to 
find the optimal solution for linear models 


Introduction: 


In this chapter, we present the modified Neutrosophic 
Simplex algorithm, which was introduced to address the 
difficulty we faced when applying the direct Simplex 
algorithm, which is the large number of calculations required 
to be performed in each step of the solution, which requires a 


lot of time and effort. 


5-1- Steps of the modified simplex neutrosophic 


algorithm: 


We explain the steps of the modified simplex algorithm 
through the following neutrosophic linear mathematical 


model: 


Max Z = Nc,x, + NcgX2 ++: + Ney Xp 


44X41 + A42X2 + ae 
24X41 + A22X2 + ae 
34X41 + A32X2 + es 


Am1X1 + Am2X%2 +°°° 


+ AynXn S Nb, 
+ AgnXn S Nbz 
+ A3znXn < Nb3 


+ AmnxXn S Nbm 


Kip sng tp a 0 


To find the optimal solution for this linear neutrosophic 
model using the modified simplex algorithm. 
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1- We write the neutrosophic linear model in standard form, 
we get the following model: 


max Z = CinX1 a ConX2 Siete 2 CunXn = Oy, ot Oy2 Heh 
+ OY, 


Ay 4Xq + Ay2X2 Fo + AynXy + V1 = Diy 
4X4 + Ag2X2 ++ + AgnXy + 2 = bon 
34X1 + A32X2 ++ + AgnXyn + 3 = ban 
Am1X1 + Am2X2 +++ AmnXn + Vm = bmn 
X41 X25 0005 Xp V1» V2» V3s +++ Vm 2 0 
2- We convert the regular linear model to the basal form and 
place the coefficients in a short table whose first column 
includes the basic variables and whose top row includes the 
non- basic variables only. 

a. We define the pivot column, which is the column 
corresponding to the largest positive value in the 
objective function row because the objective 
function is a maximization function (but if the 
objective function is a minimization function, it is 
the column corresponding to the most negative 
values). Let this column be the column of the 
variable x,. 

b. We define the pivot row; The pivot row is 
determined by the following indicator: 

bin] _ Pen 


6 = min | =—>0; a,;,>0,biny > 0 
Ais Ats 


Let this row be the line of the base variable y, 
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Then the pivot element is the element resulting from the 
intersection of the fulcrum column and the pivot row, 1.e., We 
explain the second step in the following table: 


sic variables 
x4 Xo x Nb, 
Basic Variables 
V1 A141 Gis .|\\sesene nedae | Wages Nb, 
V2 A241 Gag | xen Aon Nb, 
Vin Cia | ae. | aaetiee aaah | ipen Nby, 
NZ Ncy | NG ||| setevs saeas’ | NC | NZ NC, 


Table No. 1: Anchor element table 


We calculate the new elements corresponding to the pivot row 
and the pivot column according to the following steps: 


1. We put opposite the pivot element a,,the reciprocal of 
1 


ats 

2. We calculate the elements of the row corresponding to 
the pivot row (except the pivot row element) by dividing 
the elements of the pivot row by the anchor element a;, 

3. We calculate all the elements of the column opposite the 
fulcrum (except the fulcrum element) by dividing the 
elements of the fulcrum column by the fulcrum element 
a;; and then multiplying them by(-1) 

4. We calculate the other elements from the following 
relationships: 
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i =a a tis 
ij — Uj — Uj = 
Ats 
Qis 
Ats 
At 
Nc! = Nc; — Nc. — 
di j oa 


We get the following table: 


-_ Aijats — Atj Gis 
Nb; Q¢s NbyQis 


j — Ncjays — NesQp; 


ic Variables 
x4 X2 
Basic Variables 
y I 
V1 Qi ae. || Kaden: 
I I 
V2 A541 Go || sce 


(1) 

(2) 

(3) 
Xp Nb! 
Qin Nb, 
ay, Nb 


Table No. 2: The first stage in searching for the optimal solution 


We apply the stopping criterion of the Simplex algorithm to the 


objective function row in Table (2) below: 


Since the objective function is of the maximize type, the 
objective function row in the table must not contain any 


positive value, (but if the objective function is of the minimize 
type, the objective function row in the new table must not 
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contain any negative value), in the event that the Criterion: We 
return to step No. (3) and repeat the same steps until the 
stopping criterion is met and we obtain the desired ideal 
solution. 

We explain the above through the following example: (This 
example was presented in research [20]) 

Example: 

A company produces two types of products A, B using four raw 
materials F,, F2, F3,F,. The quantities needed from each of 
these materials to produce one unit of each of the two 
producers A,B, the available quantities of the raw materials, 
and the profit returned from one unit of both products are 
shown in the following table: 


Prodiicts Required quantity per unit Available 
quantities of the 
A “i raw materials 
Fy 2 3 [14,20] 
Fy 2 1 [10,16] 
F3 0 3 [12,18] 
F, 3 0 [15,21] 
Profit seas per [5,8] [3,6] 
unit 
Table No. 3: Issue data 
Required: 


Finding the optimal production plan that makes the company's 
profit from the producers A, B as large as possible. 


Symbolize the quantities produced from the product A with the 
symbol x,, and from the product B with the symbol xz , the 
problem will be reformulated from the neutrosophic 
perspective as follow: 
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maxZ € [5,8] X41 = [3,6] X2 
2x1 + 3x2 < [14,20] 
2x, + x2 < [10,16] 
3x. < [12,18] 
3x, < [15,21] 
x, 20,x, 20 
We apply the modified simplex algorithm: 
1- The standard form of the previous linear model is: 
max Z € [5,8] x, + [3,6] x. + Oy, + Oy2 + Oy; + Oy, 
2X4 + 3X2 ot y= [14,20] 
2X4 + X2 + y¥2 = [10,16] 
3X2 + ¥3 = [12,18] 
3x4 + V4 = [15,21] 


X4,X2 a 0, Y1, V2» V3» V4 = 0 


2- We organize the previous information in the following 
modified simplex table: 


x b; 

Basic var. . : 
Vi 3 [14,20] 
V2 1 [10,16] 


¥3 


Objective fun. 


Table No.4: Simplex table according to the modified Neutrosophic 
simplex algorithm 
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We know [a, b] < [c,d] if a < c andb < d . Therefore. 


It is clear that max([5,8], [3,6]) = [5,8] versus to the column 
of x; , meaning that the variable x, should be placed instead 
of one of the basic variables. 


Now to demonstrate which basic variables should be ejected, 
the following calculation has been done: 


[14,20] [10,16] [15,21]] [15,21] 


a a a ee 


86€Emin 


The value of @ indicates that the row versus to the variable y,, 
and the element positioned in the cross row/column is 3 which 
is the pivot element, divide the elements of the row versus to 
y4 yields. 


3- We calculate the elements of the new table using 
relationships (1), (2), (3), we obtain the following table: 


on-basic var. 
Basic var. 


y2 = 0,4 
: [0,4] 

Ys 0 [12,18] 
“1 5,7 
; [5,7] 

Object. Fun. =8 = 5 Z — [25,56] 

3’ 3 

Table No.5: Table of the first step in searching for the optimal 

solution 
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4- We apply a stopping criterion in the algorithm. We find: 
It still there is a non-negative value in the row of the objective 
function (i.e., [3,6]) 
This means that we have not yet reached the optimal solution, 
SO we repeat the previous steps as follows: 
which is versus to the xz column, this leads to the fact that the 
variable xz should be entered into the basic variables. Now, the 
question is: which basic variable should be ejected? 
Track the following calculation to answer this question, 


And here too because [a,b] < [c,d] if a<c and b<d 
» Therefore. 


oe ain en [0,4] a _ [46] 4 


a 1 * 3 


3g gral 

which is versus to the slack variable y,, the pivot element equal 
3, hence the row versus to y, should be divided by 3 , the 
required calculations yield the following tableau: 


—_ on-basic var. 


—4 
9 
¥3 : 8,12 
: [8,12] 
x4 1 [5,7] 
Object. Fun. [—6,—1] Z — [29,68] 


Table No. 6: Final solution table 
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We apply the algorithm stopping criterion. We find that the 
criterion has been met and thus we have reached the optimal 
solution. 


The optimal solution for the linear model is: 
4 4 
The value of the objective function corresponds to: 
4 
max Z € [5,8]. [5,7] + [3,6]. l= 2| = [25,56] + [4,12] = [29,68] 


It is clear from the row of the objective function that all the 
elements are neutrosophic negative numbers, this means that 
we have reached to the optimal solution is: 


4 —4 


Substitute the above optimal solution into the objective 
maximum function, yields: 

This means that the company must produce quantity xj € 
[5,7|,of product A and quantity x3 € |=, 2]of product B, 


thereby achieving a maximum profit of: 
4 
max Z € [5,8].[5,7] + [3,6]. |= 2| = [25,56] + [4,12] = [29,68] 


To compare between the modified Simplex method and the 
direct Simplex method, we solved the same example using 
the direct Simplex algorithm. Below are the solution tables: 
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Basic var. 


Vi 


V2 


3 


Objective fun. 


Table No. 7: 


X4 V1 V2 V3 Va bj 

3 1 0 0 0 [14,20] 
1 0 1 0 0 [10,16] 
3 0 0 1 0 [12,18] 


[3,6] 0 


Simplex table according to the direct neutrosophic 


simplex algorithm 


Basic var. 


on-basic var. 


—2 
3 

V3 0 0 0 0 [12,18] 
1 

X4 1 0 0 = [5,7] 
3 

: —-8 -5 
Object. Fun. 0 0 0 > ] | 2 - [25,56] 


Table No. 8: Table of the first step in searching for the optimal 


solution 
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—4 —4 
0 0 1 0 — — 

¥2 9 [ 3 , 2] 
Z 

V3 0 0 0 1 3 [8,12] 
1 

X4 1 0 0 0 = [5,7] 
é| 

, —6 
Object. Fun. 0 0 0 0 > —1] Z — [29,68] 


Table No. 9: Final solution table 


It is clear from the row of the objective function that all the 
elements are either zero or neutrosophic negative numbers, this 
means that we have reached to the optimal solution is: 


4 —A 
xj € [5,71.*3 € [5,2], 93 €|>-.2],95 € [12197 = 9% = 0 


Substitute the above optimal solution into the objective 
maximum function: 


max Z € [5,8]. [5,7] + [3.6]. 2 = [25,56] + [4,12] = [29,68] 
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Conclusion and results: 


From the previou’s study, we note that we obtained the same 
optimal solution that was obtained when we used the direct 
simplex method, but with a much smaller number of 
calculations than the number we did in the direct simplex 
algorithm, as is clear by comparing the solution tables using 
the modified simplex method, Tables No. (4)., No. (5), No. (6), 
with solution tables using the direct simplex method, Tables 
No. (7), No. (8), No. (9). Therefore, in order to shorten time 
and effort, we focus on the necessity of using the modified 
Simplex method to find the optimal solution for linear models. 
Especially when the model contains a large number of 
variables and constraints. 
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Chapter VI 


Finding a rule solution for linear models 
using artificial variables 
Introduction. 
6-1- Artificial base simplex algorithm. 
6-2- Processing the model and all constraints of type equals. 
6-3- Processing model constraints mixed. 


Conclusion. 


-85- 


Neutrosophic linear models and algorithms to find their optimal solution 


Chapter VI 


Finding a rule solution for linear models 
using artificial variables 


Introduction: 


In this chapter, we present the simplex algorithm with 
neutrosophic artificial variables, which is preferred for use in 
linear models in which there is no ready-made base through 
which we can obtain the optimal solution. Therefore, artificial 
variables are added to the constraints, the number of which is 
equal to the number of constraints that do not contain a base 
variable, and as a first step in the research. Regarding the 
optimal solution, we must get rid of all artificial variables and 
convert them to non-basic variables so that they take the value 
of zero and thus do not affect the ideal solution of the linear 
model. We explain the above through the following study: 


6-1- Artificial base simplex algorithm: 


The purpose of solving linear models is to choose the optimal 
solution from the set of acceptable solutions. This is done 
based on a base solution that is improved using the direct 
simplex algorithm, consists of three basic stages: 


1. The stage of converting the imposed model into an 
equivalent systematic form. [20] 

2. The stage of converting the regular form into a basic 
form to obtain the non-negative basic solutions. 

3. The stage of searching for the ideal solution required 
from among the non-negative basic solutions [19]. 
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Therefore, the process of searching for the optimal solution 
does not begin until after obtaining a base solution, but in 
many linear models we face great difficulty in obtaining the 
base solution, so the simplex method with an artificial base 
was proposed, where a base is formed consisting of a set of 
artificial variables that are not negativity is added to 
constraints that do not contain a basic variable, thus obtaining 
the basic solution. Then we improve it using the direct 
simplex algorithm until we obtain the optimal solution. In this 
research, we will reformulate the simplex algorithm with an 
artificial basis to find the optimal solution for linear models 
for which it is difficult to obtain a basis solution, using the 
concepts of neutrosophic. 


Text of the issue: 


Find the optimal solution for the following neutrosophic 
linear model: 


Max Z = NC, x, + NCoX2 ++: + NC) Xn + NCo 

Constraints: 

Q44X1 + Ay2QX2 te F AyynXn + & = ND, 

A24X1 + Ag2X_ $e + AgnXy + Eg = Nbp 

A34X1 + Ag2X_ $e + AgnXy + €3 = Nb3 

nay + zty +o + Oman + Em = NB 
Nipkow tye 0 
Where: 
NG =G+teE, Nb=b, +6), ay, 


PH LZ 25 1 = 125M 
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are constants having set or interval values according to the 
nature of the given problem, x; are decision variables. It is 
worthy to mention that the coefficients subscribed by the index 
Nare of neutrosophic values.The objective function coefficients 
NC,, NC3,...,NC, have neutrosophic meaning are intervals of 
possible values: That is, Nc; € [Aja Aj], where Aj1,A;2 are the 
upper and the lower bounds of the objective variables x; 
respectively, 7 = 1,2,...,n. Also, we have the values of the 
right-hand side of the inequality constraints Nb,, Nbp,..., Nbm 
are regarded as neutrosophic interval values: 


Nb; € [Hit Hiz|, here, Hj4,Ui2 are the upper and the lower 
bounds of the constraint i = 1,2,...,m. 


In the previous model, we note that the number of variables is 
n and the number of constraints is m, and this model is in the 
standard form. 


We move to the second stage, which is to find a basic solution. 
Here we use the simplex algorithm with an artificial base, 
which is represented by the following: 


1- From the standard form, we form an artificial base form 
by adding to the left side of each of the constraint 
equations a non-negative artificial variable ¢;. Thus, we 
form a base consisting of the non-negative variables 
Eo iy ey, 

2- Since the artificial variables are introduced into 
constraints that were originally linear equations, these 
variables must take the value of zero so that the linear 
constraints are not affected. 

3- Therefore, we must move all of them from the base until 
they become non-base variables, and in order to be able 
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to make this transition, we use the direct simplex 
algorithm. 

4- We introduce these variables into the objective function 
with the likes of M (where M is a sufficiently large 
positive number that is at least greater than any |N cl) 
and preceded by a minus sign (because the objective 
function is a maximization function) so that we do not 
transfer them back to the base variables again. second. 

5- We obtain the following basic form of the neutrosophic 
linear model: 

Max Z = NC,x, + NCpXx2 +++ + NC,X, — Mey — M eg —- + 
—Mé,+NCo 


Constraints: 

44X41 + Q42X2 + i + AynXn + Ey = Nb, 

24X41 + A22X2 + i + A2nXn + E> = Nb, 

34X41 + 132X2 + ais + AznXn + €3 = Nbz 

Oni hy HO Xs htt Gade Pen = Nba 
xj 20 ,€,>0,Nb; > 0;j = 1,2,..,n andi = 1,2,...,m 

6- After obtaining the basic solution, we use the direct 
simplex algorithm in order to improve this solution to 


reach the optimal solution. Therefore, we arrange the 
previous information in a table as follows: 
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ariables 
x4 Xz se Xn Ey Ey a | ea b; 

Basic 

Ey ay4 A12 ant Ain 1 0 mee 0 b, 

Ey Az Az2 he Aon 0 1 0 0 by 

Em Am1 Am2 ae Onin 0 0 as 1 bin 
Objecti 

nine NG? |G rl) | ane CRN eae tel | ads | ae 

function 


Table No. (1) General data of the model 


We get rid of the artificial variables. Here we study the 
constants b; corresponding to the artificial variables and choose 
the largest of them, let it be b, corresponding to the variable ¢; 
and we consider its row to be the pivot row. Then we 
determine the pivot element in it by dividing the elements of 
the objective function row (elements NC;) by the elements of 
the €, row and then we take the smallest positive ratio @ where: 


NC, 
Ats 


NG 
—>0 


6 = Min 
tj 


J 


Where a;; > 0, then the pivot element is a,,, and we exchange 
the variables x, and ¢,,According to the direct neutrosophic 
Simplex algorithm instructions, see [19]. We repeat step (7) 
until we get rid of all the artificial variables and obtain a 
normal base consisting of the basic variables. 


After getting rid of the artificial variables, we return to 
working according to the direct neutrosophic simplex 
algorithm. 
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6-2- Processing the model and all constraints of type 
equals: 


We explain how to find the optimal solution for linear models 
whose constraints are all equal constraints using the simplex 
algorithm with a synthetic rule through the following example: 


Example 1: 
Find the ideal solution for the following linear model: 


Max Z = —12x, + [6,9|x2 + 3x3 


Constraints: 
8x, —X2+ 4x3 = [4,6] 
6x, — 3x2 + 3x3 = [-12, —-9] 
X14,X2,X3 20 
Solution: 


1- We convert the model to the standard form, multiply the 
second equation by (-1) and we obtain the following 
model: 


Find a rule solution for the following neutrosophic linear 
model: 


Max Z = —12x, + [6,9|x2 + 3x3 
Constraints: 
8x, —X2+ 4x3 = [4,6] 
—6x, + 3x2 — 3x3 = [9,12] 


X1,X2,X3 20 
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2- We add the artificial variables and enter them into the 
objective function with a capital letter M preceded by a 
minus sign. Here we take M = 15. 


Find a rule solution for the following neutrosophic linear 
model: 
Max Z = —12x, + [6,9]x, + 3x3 —15¢, —15¢, 
Constraints: 
8x, —X2. + 4x3 = [4,6] 
—6x, + 3x, — 3x3 = [9,12] 
X1,%X2,X3,€1,€2 2 O 


We arrange the previous information in the following table: 


Variables 
x4 X3 Ey 2) bi 


Basic 


Ey 


Objective 
function 


Table No. (2) :Artificial base table 


Since the rule is artificial, we study the constants b,;and find 
that the largest of them belongs to the group [9,12| 
corresponding to the variable €,. Therefore, we divide the 
objective function row by the positive elements in the €, row 
and calculate the index @, and we find that: 


ss [6,9] 


6 = Mi 


J 
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Thus, the pivot element is (3) corresponding to x2. Therefore, 
we replace xz with €2, then the variable x, becomes a base 
variable and €, comes out of the base. We perform the 
necessary calculations and obtain the following table: 


iables 
Xs x3 Ey €2 b; 
Basic 
1 
Xo 1 —1 0 _ [3,4] 
3 
Objective 0 | [9,12] | -15 | [-18,-17] | Z— [18,36] 
function 


Table No. (3): The first change table in the base 
The artificial variable €, is still present in the base, so we 
perform another substitution, adopting the pivot line as the line 
opposite it. In order to determine the pivot column, we 
calculate the index 0, we find: 


tT [0,6] [9,12]] _ [0,6] 
6’ 3 6 

Thus, the pivot element is (6) corresponding to x,, so we move 

x, to the base instead of €,, so we get the following table: 


= Min 
j 


iables 
; xy | X2 Ey ED b; 
Basic 
se 0 |1 : 4 e 22 
3 9 Pe 
Objective 50 
function ° ° [16-19 |-13—= Z — [18,46] 


Table No. (4): The second change in the base 
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From the previous table, we note that the base variables x1, x2, 
and thus we have obtained an initial solution for the linear 
model, which It gives us the following rule solution: 


7 10 16 22 
(1 = XZ E ee X32 = 0, = 0, &2 = 0) 


But it is clear from the table that this solution is not the ideal 
solution because in the target function line there is a positive 
value corresponding to the variable x3. Therefore, we apply the 
direct simplex algorithm to improve the basic solution. We 
obtain the ideal solution from the following table: 


riables 
x4 X2 | X3 E4 2) bj 
Basic 
1 7 10 
Xg 2/o|] 1 1 z =| 
9 3° 3 
1 4 16 22 
Xp o |1] 0 es es | | 
3 9 3.3 
—53 
2480. | 0" [227224] |-19, =| Z — [39,76] 


Table No. (5): The optimal solution for the model 
Optimal solution for the linear model: 


[> Ze e[n5 10 


x, =0,x%.€ XZ E ,€, = 0,€ =0 


In this solution, the goal function takes its greatest value, which 
iS: 
Z € [39,76| 


The solution can be verified by substituting the constraints and 
the objective function statement, we note that the values in the 
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ideal solution of the previous linear model are neutrosophic 
values. 


6-3- Processing model constraints mixed: 


We illustrate how to find the optimal solution for linear models 
with mixed constraints using the simplex algorithm with a 
synthetic rule through the following example: 


Example 2: 
Find the ideal solution for the following linear model: 
MinZ = —3x, + [8,10]x, + [0,6|x3 
Constraints: 
Mp = 2Xe# xe = [3,7] 
—4x, + x2 + 2x3 = [9,6] 
2X, —-xX3=1 
X1,X2,X3 20 
Converting this model to standard form the problem becomes: 
Find the ideal solution for the following linear model: 
MinZ = —3x, + [8,10]x, + [0,6]x3 + Oy, + Oy, 
Constraints: 
Xy — 2X2 +X3+ y, = [3,7] 
—4x, + X2 + 2x3 — y, = [9,6] 
2X, —X3=1 
X4,X2,)X3,Y1,Y2 2 O 
The variable y, in the first constraint is a basic variable, and 


since there are no other basic variables, we add artificial 


-95- 


Neutrosophic linear models and algorithms to find their optimal solution 


variables to the second and third restrictions and enter them 
into the objective function in sufficiently positive times 
because the model is a minimization model, and thus we obtain 
the following basic form: 


Find the ideal solution for the following linear model: 
Min Z = —3x, + [8,10]x, + [0,6]x3 + Oy, + Oy, + 12, + 12e, 
Constraints: 
Xy — 2X2 +X + y, = [3,7] 
—4x%, +X) + 2x3 —-— yo te, = [9,6] 
2X, —X3+€=1 
X4,X2,X3, V1, V2, €1,€2 2 O 


We follow the same steps mentioned in Example 1 to remove 
the artificial variables from the base and insert the basic 
variables. After obtaining the base solution, we use the direct 
simplex method to find the optimal solution. 


Important Notes: 


1- If the row ¢; does not include a positive element and 
b, > 0, this indicates a conflict of constraints and the 
problem is unsolvable. 


, oe . NG; 
2- If we cannot find a positive ratio ae we calculate the 
tj 


largest negative ratio 0’ where: 


NG NC, 
6’ = Max |—— < 0] = 
at; ats 


Where a;; > 0, so a;, is the pivot element and it is definitely a 


positive element. 
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Conclusion: 


Through the previous study, we presented one of the important 
methods for finding the optimal solution for neutrosophic 
linear models, which is the synthetic simplex method that we 
resort to when we are unable to find a rule solution. We found 
that the optimal solution that we obtained is neutrosophic 
values, indeterminate values, perfectly defined, belonging to a 
field that represents its minimum. The smallest value that the 
objective function can take, and the highest alone represents 
the highest value of the objective function, which is 
proportional to the conditions surrounding the system’s 
operating environment, which can be represented by the linear 
model. 
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Chapter VII 


Neutrosophic Conjugate Linear Models and 
the Dual Algorithm 
Introduction. 
7-1- Neutrosophic companion models. 
7-1-1- The matrix form of the neutrosophic conjugate models. 


7-1-2- Finding neutrosophic dual models using the double 
table. 

7-1-3- Constructing neutrosophic dual linear models using 
tables. 


7-2- formulation of the binary neutrosophic algorithm. 
7-2-1- Steps of the binary simplex algorithm. 


7-2-2- Binary simplex algorithm for the original and dual 
models. 


7-3- Economic interpretation of the dual models. 


Conclusion. 
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Chapter VII 


Neutrosophic Conjugate Linear Models and 
the Dual Algorithm 


Introduction: 


In our practical life, we encounter many problems that are 
formulated in the form of linear mathematical models 
consisting of an objective function and a set of constraints in 
the form of equations or inequalities. The linear model is 
written in many formulas that are distinguished by the type of 
the objective function and the form of the constraints. The 
formulas of linear models are explained in the second chapter 
of this book, and we mentioned previously that each of these 
formulas has a use. For example, when we want to find the 
optimal solution for a linear model, we must first put it in the 
standard form. We mentioned that we use symmetric formulas 
in the dual theory, which is one of the most important theories 
in linear programming, and its basic idea is that for every 
Linear model there is a conjugate linear model, as solving the 
original linear model gives a solution to the conjugate model, 
and therefore when solving the linear programming model, we 
actually get solutions for two linear models. In this chapter, we 
present a study of the neutrosophic conjugate models and the 
dual simplex algorithm that works to find the optimal solution 
for the two models. The original and the accompanying ones at 
the same time. The importance of this algorithm is evident in 
that it is relied upon in several operations research topics such 
as, integer programming algorithms, some nonlinear 
programming algorithms, sensitivity analysis in _ linear 
programming. 
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7-1- Neutrosophic companion models: 


7-1-1- The matrix form of the neutrosophic conjugate 
models: 


To find the associated model for a given neutrosophic linear 
model using matrices, we put the neutrosophic linear model in 
the symmetrical form. As we found in the second chapter, the 
linear model is in the symmetrical form if all the variables are 
constrained to be non-negative and if all constraints are given 
in the form of inequalities (and the inequalities of the 
constraints of the maximization model must be they are written 
in the form (< less than or equal to) while the inequalities of 
the minimization model constraints must be in the form (= 
greater than or equal to), then the linear model is written in the 
neutrospheric form in one of two cases: 


The first case: The original model is symmetrical and of the 
maximization type: 


Original model: 


Find 
NZ = NC X — Max 
Constants: 
NAX<NB 
X20 

Where: 

Nay, = Nay2...NQin Nb, Ncy x4 
NaS Nay, Nag2 ...NQ2n NB = Nb, NC = Ncz Y= X2 

Naint — N@ynz Nan Nbin Nn Kn 
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The dual linear model: 


Find 


Constants: 


Where: 


Nay, 


NAT = Nay2 


Nain 


The second 


miniaturized: 


NL=NBY— Min 


NATY > NC 
Y20 
Naz,  NQmi Nb, Nc, V1 
Naz  NAm2 NB = Nb, NC = Nc» ys V2 
NaQzn - NAmn Nby, NCy Ym 


case: The model is symmetrical and 


Original model: 


Find 


Constants: 


Where: 


Nay, 


na =| N42 


Nami 


NZ =NC X— Min 


NAX = NB 
X20 
Nay we NQin Nb, Ncy x4 
Nay  NQon NB - Nb» NC = Nc» x= X2 
Naz »-» NAmn Nb», Ney Xn 


The dual linear model: 


Find 
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NL = NBY— Max 


Constants: 
NATY < NC 
Y>0 
Where: 
Nay, Naz,...Nani Nb, Nex V4 
NAT = Nay. ~=NQz2 «NA? NB = Nb» NC = Nc Y= Y2 
Nain — Naoy  NQmn Nbin Nn Yn 


We summarize the process of finding neutrosophic dual 
models using matrices in the following steps: 


1. 


We define a new non-negative variable for each 
constraint of the original model 

We make the wind (cost) vector in the original model a 
column vector of constants in the companion model 


. We make the constants column vector in the original 


model the cost (profit) vector in the companion model 
We transform a matrix of the parsimony of the variables 
of the constraints in the original model into the 
parsimony of the variables in the accompanying model 
We reverse the direction of the constraint inequalities 
We reverse the direction of the examples, that is, we 
change the increase to the maximum limit to a decrease 
to the minimum limit, and vice versa. 


7-1-2- Finding neutrosophic dual models using the double 


table: 


We previously found that we can write linear models in three 


forms: 
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The matrix form is as shown in the previous paragraph. 


The following short form: 
n 
NZ = YG zm Ej Xj — (Max or Min) 
j=1 


Constants: 


n > 
> Naijx; (=) b, +6; ; i=1,2,..,m 
j=l 


xj 20 
The detailed figure follows: 
Find 
NZ = NeyX1 + NCpxX2 +++: + NCyXy — (Max or Min) 
Constants: 
> 


NQj4X1 + NQj2X2 + +++ + NAjinXy (s) Nb, ;i=1,2,...,m 


Xs Xoo ig ty = OV 
To find the dual model, we put the neutrosophic linear model 
in the symmetrical form, and here we distinguish two cases: 


The first case: 


The original model is symmetrical and of the maximization 
type 


n 
NZ = XG + Ej xj — Max 
j=l 


Constants: 
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n 
>, Naijx; Sb 20; 20S 42.5m 
j=1 


J = 
The second case: 


The model is symmetrical and miniaturized 


n 
NL = XG + Ej )x; — Min 
j=l 


Constants: 


n 
>, Naijx; = b; + 6; F i= Ly. wey MM 


j=l 


In both cases, we have x; 20, which are the decision 
variables, unknown values that we obtain after solving the 
linear model. 
Nc; =c; + €; and Nb; = b, + 6; and Najj = aij + wij where( 
j =1,2,...,n,i=1,2,...,m) are the data of the issue under 
study, and they are neutrosophic values, unspecified values that 
enjoy a margin of freedom and are taken according to the 
nature of the situation represented by the linear model. 
7-1-3- Constructing neutrosophic dual linear models using 
tables: 

To build linear neutrosophic models using tables, we draw a 
double table for the original and utility models according to the 
following steps: 

1. The coefficients of the objective function in the original 


model are the constants column in the companion model, 
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and the constants column in the original model are the 
coefficients of the objective function in the companion 
model. 

2. We invert the signs of the inequalities of the constraints 
(if they were in the original model of type =, they 
become in the accompanying model of type =). 

3. We change the objective from maximizing in the original 
model to minimizing in the accompanying model. 

4. We place each constraint (row) in the original model 
corresponding to a column in the accompanying model, 
meaning there is one variable for each constraint in the 
original model. 

5. The variables in the original model and the dual model 
satisfy the non-negativity constraints. 


We explain the above through the following two cases: 

The first case: 

The original model is symmetrical and of the maximization 
type: 

The first case: The original model is symmetrical and of the 
maximization type 


Find 
NZ = Nex, + Ne2X2 +++» + NCyX, — Max 
Constants: 
NQ44X1 + NQ42X2 ++» +NQyynXn S ND, 
Nay4X1 + NQg2X2 +++ + NAgnXyn S Nb 


NQmyiX1 + NAm2X2 ++: + NAmnXn = Nbm 


Kent = U 
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The binary table for the original model and_ the 
accompanying model is as follows: 


Nay4X1 + Nay2X2 + °°: 
+ NaynXn 

Z NQ24X1 + NAg2X2 +°:: 

+ NagnXn 


IA 


Nami X41 + Nam2X2 + see 
+ NG Xn 


IV 


Non- X44 XQ) ey Xp 
negativity 
constraints 


Nay1y1 + Naz1y2 + ° 
+ NamiVm 
2 NaQy2V1 + NQz2V2 +°* = 
+ NQm2Vm 
= 
n NQyn¥1 + NAznV2 + °° = 
aE NaQmnVm 


Non- 
negativity 
constraints 


Table No. (1) Objective follower of the maximization type 
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The second case: The original model is symmetrical and of 
the reduction type: 
Find 


NL = Neyx, + NC2X2 +++» + NeyxX, — Min 
Constants: 
NQy4X4 + NQ42X2 + eo + NainXn > Nb, 


Naz4X4 + NQz2X2 + aa + NQznXn = Nb, 


NQmiX1 + NAm2X2 + °° + NAmnXn 2 NDm 
X74, Xo; yxy = 0 


The binary table for the original model and_ the 
accompanying model is as follows: 
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Nay4X1 + Nay2X2 +°* 
+ NQynXn 

2 NQ24X1 + NAz2X2 +°:' 

+ NAgnXy 


IV 


NQmiX1 + NAm2X2 +°°: 
+ Navy Xn 

Non- Mig Rosiay Xn 
negativity 
constraints 


IV 


Nay; + Nagi y2 ++ 
+ NamiVm 
2 Nay42V1 + Naz2V2 +°* < 
+ Nam2Vm 
< 
n Nayn¥1 + NAzn¥2 + °° < 
“Tr NaQmnVm 


Non- 
negativity 
constraints 


Table No. (2) objective follower in the original model of the reduce 
type 
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7-2- formulation of the binary neutrosophic 
algorithm. 


The binary simplex algorithm (for the original and dual 
models) neutrosophic. Through this algorithm, we can find the 
two ideal solutions for both the original and dual models at the 
same time. Before starting the binary simplex algorithm, we 
must mention the modified simplex algorithm that we will use 
within the steps of the binary algorithm. 


Modified simplex algorithm: 


In the modified Simplex algorithm, after converting the regular 
linear model to the basic form, we place the coefficients in a 
short table whose first column includes the basic variables and 
whose top row includes the non-basic variables only. We 
define the pivot column, which is the column corresponding to 
the largest positive value in the objective function row if the 
objective function is a maximization function (but if the 
objective function is a minimization function, it is the column 
corresponding to the most negative values). Let this column be 
the column of the variable x,. We define the pivot row. The 
pivot row is determined. Through the following indicator: 

Nb; Nb, 
Naj;} Nats 


N@ = min| >0; Na; >0,Nb; > 0 

Let this line be the line of the base variable y,, then the anchor 
element is the element resulting from the intersection of the 
anchor column and the anchor line, 1.e., Na,;. Then we 
calculate the new elements corresponding to the anchor line 
and the anchor column as follows: 
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1. We put opposite the pivot element Na,, the reciprocal of 
1 


Nats 

2. We calculate the elements of the row corresponding to 
the pivot row (except the pivot element) by dividing the 
elements of the pivot row by the pivot element Na;, 

3. We calculate all the elements of the column opposite the 
pivot (except the pivot element) by dividing the elements 
of the pivot column by the pivot element Na;,, and then 
multiplying them by (-1) 

4. We calculate the other elements from the following 
relationships: 


Na; Nb;Na,. — Nb,Na; 
Nb; = Nb; _ Nb, ls = l ts t us 
Nats Nats 
Na, Na;;Na,. — Na;;Na; 
Na;; = Na — Nag; 2 eee 


Nats Nats 
Na;; Nc;Na;. — Nc.Na;; 
: Nats Nats 


We apply the stopping criterion of the direct Simplex algorithm 
on the objective function row. If the objective function is of the 
maximize type, the objective function row in the table must not 
contain any positive value. (But if the objective function is of 
the minimization type, the objective function row in the new 
table must not be contains any negative value), if the criterion 
is not met, we repeat the same steps until the stopping criterion 
is met and we obtain the desired ideal solution. 


7-2-1- Steps of the binary simplex algorithm: 


a. We write the two models in basic form by adding or 
subtracting additional variables or using synthetic 
variables and isolating the non-restricting variables. 
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Basal form of the original model: 
Find 
NZ = Nex, + Ne2X2 +++ + NCynXp + Ou, + OU2 + +++ + OU, — Max 
Constans: 
NQy4X4 + NQ42X2 feet NQynXn + Uy = Nb, 
Naz4X4 + NQp2X2 Se NQznXn + Uz = Nb, 


NQm1*1 + NAm2X2 +°1° + NAmnXn + Um = Nbm 
xj 20 ;j=1,2,...,n 
u; 20 ;i=1,2,...,m 
Here we do not require that Nb; = 0. 
Basic form of the dual model: 
Find 
NL = Nbyy; + Nbogy2 +++: + NDiVm + OV, + OV2 +++» + OV, — Min 
Constans: 
NQy1¥1 + Naz, y2 + °° + NAmi1 Ym — V1 = Ney 
Nay2¥1 + NA22Y2 + °° + NAm2Vm — V2 = Nez 


Nayn¥1 + NAgn¥2 + ++ + NAmnYm — Vn = Ney 
y, 20 5i=1,2,..,m 
vy, 20 5j =1,2,...,n 
Here we do not require that Nc; 2 0 
The two models have the same coefficients, and the matrix of 
instances of the dual model is the transpose of the matrix of 
instances of the original model. We write the two models in the 


following binary table: 
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NQ44X1 + NQ42X2 +++: + NQyyXy 
+ uy, 

2 NQp4X1 + NAz2X2 +1 + NQgnXy | = 

+ uz 


NQy1X1 + NAm2%X2 +°°° 
NG ae ie 

Non- XX oq, ey Xy Uy, Ug, +, Um 2 
negativity 
constraints 


Nay1y1 + Naziy2 + °° 

+ NQmi¥m Vy 
2 Nay2y1 + Nazzy2 + = 
a NQm2Ym — V2 


n NQyn¥1 + NAgn¥2 ++ = 
+ NQmnYm Vy 


Non- 
negativity 
constraints 


Table No. (3) Standard format for the original and companion 
models 
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b. We place the variables and coefficients of the original 
model in the modified simplex table, and we place the 
variables of the dual model outside the table as follows: 


Table No. (4): The binary table for the original and dual models 
according to the modified Simplex algorithm 
7-2-2- Binary simplex algorithm for the original and dual 
models: 


From the modified simplex algorithm of the original model, we 
obtain the optimal solution of the original model when all the 
elements are in the last row (the objective function row of the 
original model) Nc; < 0 ;j = 1,2,...,n and at the same time 
all the elements are in the column The last (associated 
objective function column) Nb; = 0 ;i = 1,2,...,m and we get 
the optimal solution for the dual model when all elements in 
the last column (associated objective function column) are 
Nb; = 0 ;i =1,2,...,m and at the same time the last row (the 
objective function row of the original model) Nc; <0 ;j = 
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1,2,...,n (because it will correspond to Nc; = —v;) which are 
the two conditions Same for both models. Therefore, when 
searching for the optimal solution for both models together, we 
must work to make all elements Nb; => 0 ;i = 1,2,...,m and to 
make all elements Nc; < 0 ;j = 1,2,...,n to achieve this, we 
rely on one of the two models, put its variables and coefficients 
in a table, and place the dual model in an external frame of that 
table. In general, we find that the necessity of placing the two 
models in a short table does not allow us to get rid of the 
negative constants on the right side, and therefore the general 
case of the previous binary table can it must include negative 
constants Nb; <0 , and the elements of the last row can 
include positive elements Nc; > 0 , so when searching for the 
optimal solution for the two models, we must work to address 
these elements based on one of the two models. 


Depending on the original model, we do this in two stages: 
The first stage: 


We make the constant Nb; non-negative, which corresponds to 
obtaining a non-negative basic solution for the original model. 


The second stage: 


We make all elements of the objective function row non- 
positive (in the case of the objective function, maximization), 
and this corresponds to obtaining the optimal solution required 
for the original model. 


Based on the dual model, we do this in two stages: 
The first stage: 


We must make the elements of the dual model objective 
function column 


eee 
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Nb; = 0 ;i=1,2,...,m_ ,The last row is non-negative 
The second stage: 


We must make the free constants for the dual model 
—Nc; non-positive, and this corresponds to obtaining the 
optimal solution for the dual model. We explain the above 
through the following example: 


Find the optimal solution for both the following neutrosophic 
linear model and its dual using the binary algorithm 


Example1: 


Find: 
[5,8] x, + [3,6] x. — Max 


Constans: 
2x, + 3x2 < [14,20] 
2X, +x, < [10,16] 


3x. < [12,18] 
3x, < [15,21] 


x, 20,x, 20 


We form the binary table of the model and the dual model: 
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negativity 
constraints 


2y1 + 2y2 + 3y4 


3Y1 + ¥2 + 394 


Non- 
negativity 
constraints 


Table No. (5) The original model and its dual mode 
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In the following table, we wrote the two models in standard form: 


objective 
function 


| 2X, + 3X2 + Uy 

2 2x, +x2 +U2 = 

3 3x2 + U3 = 

ae 3x, + uy, = 
Non- X4,X2, Uy, U2, Ug, U4 > 


negativity 
constraints 


objective 


2y1 + 2y2 + 34 — V1 


3Y1 + ¥2 + 3BY4 — V2 


Non- 
negativity 
constraints 


Table No. (6) Standard format for the original model and the dual model 
We notice from the table that the standard form of the original 
model includes a ready-made base of additional variables 
U4, Uz, U3, U4, but for the dual model there is no ready-made 
base. Therefore, we multiply the two restrictions by (-1) and 
we obtain the basic form of the dual model. 
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The following table shows the basic form of the original and 
dual models: 


objective 
function 


1 2X, + 3x2 4+ Uy 
2 2X4 + X2 + uz = 
3 3X2 + U3 = 
4 3x4 + U4 = 
Non- X4,X2, U4, U2, U3, U4 = 
negativity 
constraints 


objective 
function 


—2y1 — 2y2 — 3Y4 + Vy 


—3Y1 — V2 — 394 + V2 


Non- 
negativity 


constraints 


Table No. (7): The basic shape of the original model and the dual model 


We put the two models in the modified Simplex algorithm 
table and we get the following table: 


- 118 - 


Neutrosophic linear models and algorithms to find their optimal solution 


1 


y2 


¥3 


V4 


Non-basic vibrable Dual model 


Non-basic objective 
vibrable ee 
basic vibrable 
NB; 
2 3 [14,20] 
2 1 [10,16] 
0 3 [12,18] 
3 0 [15,21] 
L-0O 
u u u 
1 2 3 U4 
Non-basic 
vibrable 
basic vibrable 
V1 —2 —2 0 —3 
V2 3 —1 —3 0 
objective [14,20] | [10,16] | [12,18] | [15,21] n= 
function 
Dual model L—0 
NB; 


Table No. (8): The binary table for the original and dual models according 
to the modified Simplex algorithm 
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The first stage: 
1- For the original model: 

Since the values in the constant’s column are all positive, we 
study the values in the objective function row and determine 
the largest positive value. We find: 

max([5,8], [3,6]) = [5,8] 
It is an expression of the variable x,. This means that it will 
enter the base. To determine the element that will exit from the 
base, we calculate the index N@, where: 

[14,20] [10,16] [15,21]} [15,21] 

2 De Be 

We find that the pivot column is the column of the non-base 


N@ Emin 


= [5,7] 


variable x,, meaning that the variable x, will enter the base 
instead of the variable u,, and the pivot element is the element 
resulting from the intersection of the pivot row and the pivot 
column, which is (3) 
We perform the switching between variables using a modified 
simplex algorithm. 

2- For the dual model: 

We study the elements of the objective function row. We 
notice that all the values are positive. Therefore, we study the 
elements of the constant’s column. We find that they are all 
negative values. We choose the most negative of them, which 
is (—[5,8]) which is the row of the base variable v, so its row 
is the pivot row. To determine the pivot column and the pivot 
element, we calculate the index N@’where: 

[14,20] [10,16] [15,21] [15,21] 


N@’ € Max 257 ag es 3 


So, the column of the non-base variable u, is the pivot column, 
meaning that the variable u, will enter the base instead of the 
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variable v,, and the pivot element is the element resulting from 
the intersection of the pivot row and the pivot column, which is 
(-3). We perform the switching between the variables using the 
modified simplex algorithm, from (1) and (2) We get the 
following double table: 


V4 V2 
Non-basic objective function 
vibrable Original model 
NB; 
basic vibrable 
¥1 is 3 [4,6] 
3 
o 
ra = 
£ |v = 1 [0,4] 
' 3 
2 
& | ¥3 0 3 [12,18] 
5 
Zz 
Vy 1 0 [5,7] 
= [25,56] 
Uy U2 U3 x4 
Non-basic 
vibrable basic 
able 
Va Z 2 0 as 
3 3 3 
V2 —3 —1 =5 0 
objective function [4,6] [0, 4] [1 2; 18] [5,7] 
Original model a [2 
Ne i : 


Table No. (9): The binary table for the first stage, the solution 
according to the original and dual models 
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The second phase: 
We apply the stopping criterion of the algorithm 
For the original model: 
Since the values in the constant’s column are all positive, we 
study the values in the objective function row. We notice that 
there is a positive value, which is [3,6], meaning that we have 
not yet reached the optimal solution. Therefore, we specify the 
pivot column, which is the column of the variable x, 
corresponding to the only positive value in the objective 
function row. [3,6] In order to determine the pivot row and the 
pivot element, we calculate the index N@, where: 

[4,6] [0,4] [12,18] [4,6] 


N@ Emin “eae 8 a 


It corresponds to the base element u,, so its row is the pivot 
row and the pivot element is (3). We swap between the 
variables using the modified simplex algorithm. 

For the dual model: 

We study the elements of the objective function row. We 
notice that all the values are positive. Therefore, we study the 
elements of the constant’s column. We find that there is a 


single negative value, which is (—[3,6]), which is the line of 
the base variable 172, so its row is the pivot row. To determine 
the pivot column and the pivot element, we calculate the index 
N@' where: 
[4,6] [0,4] [12,18]] [4,6] 
N@' € Max |——,— —. 
leg oda ee =3 


So, the column of the non-base variable y, is the pivot column, 


meaning that the variable y, will enter the base instead of the 
variable vz, and the pivot element is the element resulting from 
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the intersection of the pivot row and the pivot column, which is 
(—3). We perform the switching between the variables using 
the modified simplex algorithm, from (1) and (2) We get the 
following double table: 


—Y4 Yi 
on-basic objective 
vibrable function 
Dual model 
basic vibrable NB; 
V2 — i E 2| 
4 9 3 3’ 
= 
E | y 2 = H 2| 
ia 9 3 3 
a | ¥ a =a [8,12] 
B 3 
aT 2) 1 0 [5,7] 
— [29,68] 
x2 U2 Ug x4 
on-basic 
vibrable 
basic Vibrable 
V4, 2 4 —2 —1 
9 9 
oA a1 1 1 0 
3 3 
objective E 2 E 2 [8,12] [5,7] 
function Q? 2? = 
Dual model 2 3 E 68] 
NB; 


Table No. (10): The binary algorithm table for the second stage 
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We apply the stopping criterion of the algorithm: 


1- For the original model, we study the elements of the 
objective function row until the criterion for stopping the 
algorithm is met, which is the absence of any positive 
element 

2- For the dual model, we also study the elements of the 

constants column until the criterion for stopping the 

algorithm is met, which is the absence of any negative 
element 

We find that the criterion has been met and thus we have 

reached the optimal solution 


x 


The optimal solution of the original model is: 
4 4 
x5 E 5-2 | ,u5€ |= -2| ,u3 € [8,12] ,xj € [5,7] ,uj =u, =0 


The value of the objective function corresponds to: 
NZ* = Max NZ € [29,68] 
The optimal solution of the dual model is: 
yi € [1,2] ,y4¢ €[L6] ,02 =y2 =y3 =] = 0 
The value of the objective function corresponds to: 
NL* = MinNL € [29,166| 
We note that: 
NZ* = Max NZ € [29,68] < NL* = MinNL € [29,166] 


This solution is acceptable according to the following 
theory: 


If (%1, Xz, ..., Xp,)is an acceptable solution to the original model 

of type Max and 

(V1, ¥, ++» ¥m) was an acceptable solution for the dual model of 
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type Min, so the value of the objective function of the original 
model does not exceed the value of the objective function of 
the dual model for these two solutions, that is, it is 


n m 
jJ=1 i=1 


This is for all acceptable solutions for both models (including 
the optimal solution). 


7-3- Economic interpretation of the dual models: 


We illustrate the economic interpretation of the dual model 
through the following example: 


Example2: 


A factory wants to transfer its products from two warehouses to 
three retail centres at the lowest possible cost. The following 
table shows the data provided by the factory official: 


By By B; Available 
Stores quantities 
A; [1,3] | [2,4] | [0,3] 300 
Ay [4,6] | [1.4] | [1,5] 600 
Quantities 200 300 400 900 
required ane 


The plant manager's request was for a transportation plan with 
a minimum cost so that the distribution centres’ orders could be 
met from the available quantities 


The previous issue is a balanced transfer issue because 
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2 3 
D a; = > bj = 900 
=1 j=1 


To formulate the Sthedlatieal model 
We assume x;; the quantity transported from store i where 


i = 1,2, to distribution center j, where j = 1,2,3. Thus, we 
obtain the following linear model: 
Find 
L € [1,3]x41 + [2,4]x42 + [0,3]x13 + [4,6]x2, + [1,4] x22 
+ [1,5]x.3 — Min 
Constans: 
X44 + X42 + X13 < 300 
X44 + X21 + X23 < 600 
X44 + X21 = 200 
X12 + X22 = 300 
X13 + X23 = 400 
xj20; i=1,2 ,j=1,2,3 
We write the model in the following symmetrical form: 


Since the objective function is a minimization function, all 
constraints must be of type greater than or equal to, so the 
model is written in the following symmetric form: 


Find 


L € [1,3]x41 + [2,4]x12 + [0,3]x13 + [4,6]x21 + [1,4]x22 
+ [1,5]x.3 — Min 


Constans: 
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X41 — X12 — X13 = —300 
X44 — X21 — X23 2 —600 
X41 + X21 = 200 
X12 + X22 = 300 
X13 + X23 = 400 
xj20; i=1,2 ,j=1,2,3 
Forming the dual model, we obtain the following linear model: 
Find 
Z = —300y, — 600y,z + 200y3 + 300y, + 400y, — Max 


Constans: 


V1» V2» V3» Va V5 = 0 


We formulate an appropriate text for the accompanying model 
based on the text of the original problem: 


It is clear from the original model that the factory's goal is to 
transport all of its products at the lowest possible cost: 


Text of the issue dual to the attached form: 


A transport company submitted to a factory an offer that it 
would transport the entire quantity in the first warehouse, 1.e., 
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300 units, at a price of y; monetary unit per unit, and transfer 
the entire quantity available in the second warehouse, 600, at a 
price of yz monetary units per unit. The company pledged that 
it would deliver (200,300,400) units. To the three retail 
centers, respectively. These units are sold in these centers at a 
price of (V3, ¥4, Ys) monetary units, respectively. So that you 
can convince the official in the factory that if he accepts her 
offer, the cost of transportation in his factory will be less than 
the cost if he carries out the transportation process, so that he 
carries out the transportation process. 


She used the constraints in the dual model and conducted the 
following discussion: 


You pay the cost of transporting one unit from the first factory 
to the first sales center, an amount whose value belongs to the 
range [1,3], but if you use the transport company, the cost is 
(v3 —¥1), and we have from the first entry in the 
accompanying model 


y3—-¥1 S [1,3] 
Here the official in the laboratory will notice that the 
transportation company’s offer is an appropriate offer. 


In the same way we discuss all the limitations of the dual 
model, the conclusion that the factory official will reach is that 
the cost of transportation on any route if the transportation 
company’s offer is accepted is less than or equal to the cost that 
he would pay if he himself carried out the transportation 
process. 


The transport company will adopt the values (v1, V2, V3, V4, Vs) 
,because it will achieve maximum profit through them, as the 
transport company’s profit is calculated from the relationship: 
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—300y, — 600y, + 200y3 + 300y, + 400y. 


It is the same as the objective function of the dual model, 
meaning that the dual model represents the transportation 
company that is trying to maximize its profits 


The basic theorem of association states that the optimal values 
of the model and the dual model are always equal. The 
manufacturer does not save any money because he will pay the 
transportation company the minimum cost of transportation, 
but it saves the trouble of solving the original model to 
determine the minimum cost of transportation, and for the 
transportation company, it has guaranteed the deal to transport 
the goods with the maximum profit. 


Conclusion: 


From the previous study, we arrived at a solution for the 
original and utility models at the same time, which are 
neutrosophic values from which we know the minimum and 
maximum profit that we can obtain, because the interpretation 
of the optimal solution for the original model is that it gives us 
the best production plan that makes the value of that production 
as large as possible, within Available capabilities. As for the 
optimal solution for the dual model, it gives us the best values 
for the prices of raw materials, which, if used without waste, 
will also give us the best production plan, and the result is the 
maximum profit. 
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Chapter VIII 


Some applications to neutrosophic linear 
models 
Introduction. 
8-1- Problem of the composition of mixtures. 
8-2- Problem of product mixture. 


Conclusion. 
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Chapter VIII 


Some applications to neutrosophic linear 
models 


Introduction: 


The linear programming method is one of the most used 
operations research methods in most practical fields. This 
method depends on converting the issue under study into a 
linear mathematical model, and then we find the optimal 
solution using special algorithms to solve linear models. This 
solution helps the decision makers responsible for managing 
the system. Which works according to this model to make 
sound decisions based on scientific foundations. The most 
important stage in linear programming is the stage of creating 
the linear model, that is, expressing the situation under study in 
mathematical relationships. In order to formulate the linear 
model, the following basic elements must be present: 


1. Determine the goal quantitatively, and it is expressed by 
the goal function, which is the function for which the 
maximum or minimum value is required. That is, we 
must be able to express the goal quantitatively, such as if 
the goal is to achieve the greatest profit or achieve the 
lowest cost. 

2. Determine the constraints: The constraints that express 
the available resources must be specific, that is, the 
resources must be measurable, and expressed in a 
mathematical formula in the form of inequalities or 
equals. 

3. Identifying the different alternatives: This element 
indicates that the problem should have more than one 
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solution so that linear programming can be applied, 
because if the problem had one solution, there would be 
no need to use linear programming, whose benefit is 
focused on helping to choose the best solution from 
among the acceptable solutions. 


In this chapter, we present some problems that lead to 
neutrosophic linear models, that is, we will take some or all of 
the problem data as neutrosophic values. 


8-1- Problem of the composition of mixtures: 


By mixtures, we mean anything that is installed from a number 
of materials such as diets - medicine - any metal mixture- and 
here the stretch loop is to choose the materials that enter the 
composition of this mixture so that the cost of production is as 
little as possible, the goal of putting forward this model in the 
field of education is that the student can link between 
neutrosophic equations and linear inequations as well as the 
neutrosophic function and problems from real life. 
General text of the problem: 
We want to install a mixture of raw materials Aj, Ag, ..., A,and 
the price of one unit of each of them is equal to 
NC,,NC,,...,NC, respectively, and the meal must include an 
amount of important elements 
B,,Bz,...,Bm that the quantity of each element shall not be 
less thanB,, Bo, ..., Bm 
Nb,,Nbz,..,Nb,, Unit in the order required to find the 
necessary amounts of each of the materials 
A,,A2,...,4, Which must be included in the mixture so that 
its cost is as low as possible, knowing that the content of each 
of the materials of each of the elements is shown in the 
following table: A,,Az,..,Ay, , By, Bz,..-,Bm 
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Ay ~ ||| saesdaee Ap Ay w materials 
Elements 
Nee a2 a41 By 
Nb». oo A22 a21 By 
Mia Nese, eee Am2 Qm1 Bm 
m 


Ifand Ncoj =cj te j = 1,2,...,n where is indefinite and can 


be Ore; 


& = [Aj Adil 


EE {Ayj A2;} or... else 
Also values that express the quantities of elements that must be 
available in the mixture and 
Nb, = b; +6; it = 1,2, ...,m where it is indefinite and can be 
or or ... otherwise 6; 6; = [esi Hail 6; € {eri Hai}- 
Building the Mathematical Model: 
symbolize the required amounts 


We 


materialsA,, Ao, ... 


An 


BX, Xo, ++ 


Xn 


information in the following table: 


of each of the 
and put all the 


Minimum ri A, Ay materials 
Ingredients - Elements 
o---n22---------- a a 
Nb, Be) 12 it B, 
w 22-22-2222 ----- a a 
Nb, ee 22 21 Bs 
won nn n--nn------ a a 
Nby, Amn m2 m1 Bn 
| (eee Nc» Nex Price per unit 
Xp | recor eee eeeeeeees X2 X41 Required 
Amounts 
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What is required in the problem is to determine a value for 
each of the variablesx,,X2,...,X, so that the target function 
takes the smallest value, within the conditions imposed 

Based on the data of the problem, the objective function is 
written in the following form: 


NL = NeyxX1 + NC2X2 ++: + NCpXy 


Mathematically formulate the terms we provide the following 
clarification: 


Each unit of the material A,gives us a,, unit of the 
elementB, , and thus we find that x, unit gives us a,,xX,unit 
of the element B,, and so we find for the rest of the materials 
and elements, and therefore the condition related to the element 
B,is as follows: 


44X41 + A42X2 + ae + AynXn = Nb, 


We follow in the same way for all materials and all elements 
we get the following mathematical model: 


MinNL = Neyx1 + NceoX2 ++: + NCyXy 
Within the conditions 
44X41 + Q42X2 + rata + AynXn => Nb, 


A24X4 + A22X2 + aes + A2nXn = Nb, 
Gig Xt Am Xo b+ QmyX, 2 Nb 


Kiegiagh pe OU 


It is written in the following abbreviated form: 


Ha 
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n 
Min NL = »G + &) x; 


j=1 
Subject to 
ya, *, 2b t0.- a= 12.0 
=I 
x, 20 
Where.cj te; , Db £6; ay, J = 1,2,..,n,t = 1,2,..,mx; 


are constants having set or interval values according to the 
nature of the given problem, are decision variables 


Example: 


A school wants to provide breakfast for students consisting 
four types of food that are A,, A ,A3,A, and that the price of 
one unit of each of them respectively Nc,,Nc2,Nc3,Nc, , and 
suppose that it is required that the meal includes a certain 
amount of important nutrients: —proteinsB,starch B, 
carbohydrates B3 that the amount of protein in it is not less 
than Nb, 


unit and the amount of carbohydrates is Nb, unit and the 
amount of carbohydrates is Nb3 unit, it is required to find the 
necessary amounts of substances that must be included into the 
meal so that its cost is as low as possible and contains at least 
the minimum number of nutrients required to be provided in 
the meal, if you know what each of the materials and the 
important nutrients they contain are required are shown in the 
following table: 


where N GHG + Ej and j = 1,2,3,4 where gj is indefinite 


and canbe &€ = [Ary Ap; lor 
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Ee {Ayj A2;} or... else 


Also values that express the amounts of nutrients that must be 
and i = 1,2,3 where 6; it 
is indefinite and can be 6; = [oa ‘Hil or 6; € {eri Hai} Or... 


available in the mealNb; = b; + 6; 


otherwise. 


We denote the required amounts of each of the materials 
A,,Az,A3,A, with symbols x1, %X2,X3,X4 respectively put the 
information contained in the text of the problem’s table as 


follows: 
Minimum Ay A; A, Mi “ 
Ingredients a 
Nb, Ara a13 a2 a1 By 
Nb, Ce Q23 Q22 Q21 B, 
Nb; Asa 233 432 431 Bs 
Neg Nc Nc» Ne, Price per 
unit 
X4 X3 Xy X4 Required 
Amounts 
Follow the goal: 


From the previous table we find 


Nutrient conditions: 


Nutrient protein requirement By, 


Requirement of starch nutrient B, 


Requirement of carbohydrate nutrient B 


NL = Nex, + Ne2X2 + Ne3X3 + NC4Xq 


Q44X1 + Ay2%X%2 + Ay3X3 + Ay4X, = Nb, 


24X14 + Ag2X2 + Az3X3 + Ag4X4 = Nb 
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A34X1 + A32X2 + A33%X3 + Az4X, = Nb3 
Non-negativity condition: 
X1,%X2,X3,X4, = 0 
That's when the right mathematical model becomes. 
Find the smallest value of the function 
MinNL = Nceyx, + Nc2X2 + Ne3x3 + NC4XxX4 
Within the conditions 
Q44X1 + Ay2%X%2 + Ay3X3 + Ay4X, = Nb, 
C24 X14 + Ag2X2 + Az3X3 + Ag4X, = Nb 
34 X1 + Ag2X2 + Ag3X3 + Ag4X4 = Nb3 
X4,X2,X3,X, =O 


To obtain the optimal solution, we use the neutrosophic 
simplex method described in the research 


8-2- Problem of product mixture: 


What is meant by the mixture of products is the production of 
products that include in the composition of a number of raw 
materials in all production institutions in order for the 
workflow to be ideal and achieve the maximum profit must be 
based on a scientific study through which the quantities that 
must be produced from each product of the products are 
determined by using the available resources ideally so that the 
market need is secured and a profit is made. This model can be 
presented to students as an example of the use of linear models 
then students realize that pens and notebooks and benches, 
tables, transportation-------- and other products that they use in 
their daily lives have been manufactured using process research 
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methods that depend on building a mathematical model and the 
optimal solution for this model is the ideal plan that the 
institution must adopt. 


General text of the problem: 


A production institution that can produce products 

A,, Az, ..., Ay and includes in its composition of raw materials, 
the quantities used of each B,, Bp, ..., By, of the raw materials 
in each of the products are shown in the following table: 


ae eee A> A, aterials 

Elements 

a | Kea cea a12 Qy4 B 1 
1n 

PAI ica a a22 a21 By 
2n 

ie, Wier Neen ee Am2 am1 Bm 
mn 


The quantities available to the institution of these raw materials 
are Nb,, Nbz,...,Nb,, in order where 


Nb; = b; + 6; And where it is indefinite and can be or or ---- 
otherwise, iti = 1,2,...,md; 


Oj = [Mai Hail 6; € {Hai Hoi} is required to find the amount of 
what must be produced from each of the products, knowing 
that the profit returned from one unit of each of the products is 
respectively ---- j = 1,2,..,n&€ where and 


Ay, Ap, «.,AnNCy, NC, ..,NC, Nc = cj + & 
where it is indeterminate and can be 
Ej = [Arj Aaj ee {Ayj A2;}- or other 
Building the Mathematical Model: 
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We code the quantities produced from each of the 
productsA,,A2,...,4, Bx,X2,..,X, and put all the 
information in the following table: 


Available Ae Wome thane, Az Ay 
Quantities ‘ 
wenn eee enenee a a 
Nh; ts 12 11 
wenn eee cnneee a a 
Nb, Aan - - B, 
wenn eee cnneee a a 
Nby Amn ms ee Bm 
Wes. ' |W asrecesasss Neo Nc, Profit per 
poe unit 
ee, |b Mbeeteeaanend Xo Xy Quantities 
= produced 


What is required in the problem is to determine a value for 
each of the variablesx,,X2,...,X, so that the target function 
takes the greatest value, within the imposed conditions. 

Based on the data of the problem, the objective function is 
written in the following form: 


NZ = NeqX, + NC2X2q ++: + NCyXy 


Mathematically formulate the terms we provide the following 
clarification: 


To produce one unit of the product A, , we need a,, unit of 
the material B,, and thus we find that x, unit of the productA, 
needs A,a,, unit of the materialB,, and so we find for the rest 
of the products and materials, and therefore the condition 
related to the material is as follows:B,x,A,Q11x,B,B, 
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44X41 + A42X2 + Pe + AynXn < Nb, 


We follow in the same way for all products and all materials 
we get the following mathematical model: 


Find the maximum value of the function 
MaxNZ = Nex, + Nce2X2 ++: + NCnXy 
Within the conditions 
Q44X1 + Ay2QX2 t+ AyynXn S Nd, 


24X41 + A22X2 + Re + AgnXn < Nb, 


pik TiGyaks br Gynt, = ND, 
Moy gy 


It shall be written in the following abbreviated form: 


n 
MaxNZ = XC + &)) x; 
j=l 


Subject to 


ae 
& 
~. 
II 
as 
N 
wW 


n 
> Aj jx; < b; 
j=l 


Where Cj a &j ; b; + 6;, Ajj, yi = A wy Nyl = LZ, wey MM 
are constants having set or interval values according to the 
nature of the given problem, x; are decision variables 
Example: 

A factory for manufacturing pens produces four types 


Sys Se 5Soue and uses the following raw _ materials 
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M,,M,,M,for this. The factory management wants to study 


the optimal organization of production during a period of time 
(for example, a month) and determine the monthly production 
of each product in order to achieve a maximum profit, knowing 
that the profit is directly proportional to the number of units 
sold of products. We explain the available quantities of raw 
materials required for each product and the profit returned in 


the following table: 
Thequantities 5; S3 52 51 a 
available 

Nm, is a43 a42 ay4 M, 

Nm, Aa a23 a22 a21 M, 

Nm, a4 433 a32 a31 M; 


To construct the mathematical model , we assume x ; counting 


units produced from S, 


X2 number of units produced fromS, 


X3 number of units produced fromS3 


xX, Number of units produced from S, 


During the productive period (for example, a month) we put 


the information in the following table: 


stuffs 
Available S4 S3 52 “ of 
si Nutri 
Quantities 
a a a 
Nm, Aras 13 12 11 M, 
a a a 
Nm, doa 23 22 21 M, 
a a a 
Nm 34 33 32 31 M; 
Ncg Nc3 Nc Nc, Profit per unit 
X4 X3 Xo X4 Ingredients 
produced 
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From the table we can see that the primary material condition 
M;: 

4 4X1 + Ay2X2 + A43X3 + Aygx, S Nm, 
Initial material requirement M>: 

24X14 + Az2X2 + A23X3 + Ag4X, S NM, 
Initial material requirement M;;: 

34X1 + Az2X2 + A33X3 + Az4X, < Nm 
In addition, the quantities produced must be non-negative, 1.e.: 

X4,X2,X3,X, =O 

We now define the target function If units of the same type are 
produced x, ,%3,X,,X, respectively, the profit during the 
production period will be: 

NZ = NeyXxX1 + Ne2X2 + NC3X3 + NC4X4 
Therefore, the mathematical model of the problem 1s: 
Find the maximum value of the function 


MaxNZ = Ncyx, + Nc2X2 + Nce3x3 + NC, 

Within the conditions 

44X41 + Ay2X2 +. A43X3 + Ay4X4 SNM, 

Ap4X1 + Ag2X2 + Az3X3 + Ag4X4 S< Nm, 

34X41 + Az2X2 + A33X3 + Az4X4 < Nmz3 

X4,X2,X3,X, =O 

We will apply the above to the model of optimal use of 
agricultural land, using the concepts of neutrosophic 
science. We will take data that is affected by the 


surrounding conditions, neutrosophic values. 
Text of the issue: 


aye 
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Let us assume that we have n agricultural areas (plain or 
cultivated), the area of each of which is equal to 
Ay, Az," ** ,A, , We want to plant it with m types of 
agricultural crops to secure the community’s requirements for 
it. Knowing that we need of crop i the amount b; , if the 
average productivity of one area in plain j of crop i is equal to 
Na,jtons/ha. Where j = 1,2,...... ,n andi = 1,2,..... ,m, and 
the profit returned from one unit of crop i equal to Np;, 
WhereNp; is a neutrosophic value, an undefined non-specific 
value that designates a perfect and can be any neighbor of the 


value a;;, also Np;which can be any neighbor of p;. 


p 
Required: 

Determine the amount of area needed to be cultivated with 
each crop and in all regions to achieve the greatest possible 
profit and meet the needs of society. 

Formulation of the mathematical model: 


We symbolize by x;; the amount of area in area j that must be 


ij 
cultivated with crop , and we place the data for the problem in 


the following table: 


Regions 1 rn n Order amount | profit amount 
Crops b; Npi 
1 Nay, Nay2 | 0 Nain by NP 
X11 X12 Xin 
2 Naz, IDS «|||; eh Nazn b2 NP2 
X21 X22 X2n 
m NQmi NQm2 0 | Namn Din NPm 
Xm1 Xm2 Xmn 
Available ay Az | ocr An 
space 
Qj 


Table No. (1) Issue data 
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Then we find that the conditions imposed on the variables x;; 


are: 
1- Space restrictions: 


The total area allocated to various crops in area J must be equal 
to a; , that is, it must be: 


X41 + X12 foireeeee + Xm1 = ay 
X12 + X22 + eee cee + Xm2 = az 
Xin + Xon + baie eae + Xmn = An 


2- Conditions for meeting community requirements: 


The total production of crop i in all regions must not be less 
than the amountb; , that is, it must be: 


NQ44%14 + NQ42X42 + i ea + NQynX1n = b, 
NQp4X241 + NQpz2X22 + pa + NQpznX2n = bs 
NQmixXm1 + NAm2Xm2 t-°°" Pt NGaiX an SD 


Find the objective function: 


We note that the profit resulting from the production of crop i 
only and from all regions is equal to the product of the profit 
times the quantity and i.e: 


Np (Nay, Xj, + Naj2Xj2 +++ + NainXin) 


Thus, we find that the objective function, which expresses the 
total profit resulting from all crops, is equal to: 
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From the above we get the following mathematical model: 


Find the maximum value of 


Xqq $F XyQ Fore + Xm1 = ay 
X12 + X22 + sere neres + Xm2 = Az 
Xin + Xon + Brag tape + Xmn — An 
44X41 + Ay2X%y2 Tere + Ayn Xin 2 by 
24X21 + Ag2Xo2 Fee + AgnX2n = bz 
Qm1ix*m1 ah QAm2*m2 a i + QAmn*mn = bm 
xij 20 of eae ea? eee ,m,j = 1,2, wtih yn 


Example: 

Let us assume that we want to exploit four agricultural areas 
A, Az,A3,A4, and the area of each of them, respectively, 
is 60,150,20,10, by planting them with the following crops: 
wheat, barley, cotton, tobacco, and beet, from which we need 
the following: 800, 200,600,1000,2500 Let us assume that the 
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regions’ productivity of these crops and their prices are given 
in the following table: 


egions | Aj, A Az A, | the order | Price per ton 
Crops 
wheat {4,6} 4 3 6 2500 {1400,1600} 
barley 7 5 4 {3,5} 1000 {900,1100} 
cotton 4 {9,11 } 8 5 600 {4500,6000} 
tobacco 6 {2,4} 0 0 200 {4000,5000} 
beet 3 {10,14} 10 6 800 {400,700} 
Space 60 150 20 10 
Table No. (2) Example data 
Required: 


Formulate the mathematical model for this issue so that the 
production value is as large as possible 
To formulate the mathematical model, we extract the following 


linear conditions: 


Space restrictions: 


X41 + X24 + X34 + X41 + X54 = 60 


X12 + X22 + X32 + X42 + Xs2 = 150 


X13 + X23 + X33 + X43 + X53 => 20 


X44 + X24 + X34 + X44 + X54 => 10 


Order restrictions: 
{4,6}x11 + 4x12 + 3x13 + 6x14 = 2500 
7X21 + 5X22 + 4X03 + {3,5}x24 = 1000 
4X3, + {9,11 }x32 + 8x33 + 5x34 = 600 


6X41 + {2,4}x42 + 0x43 + 0x44 = 200 


3x5, + {10,14}x.5 + 10x53 + 6x54 = 800 
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Non-Negative restrictions: 
xij 20 5t=1,2,3,4,5 and j = 1,2,3,4 
Objective function that expresses the value of production 
is: 
Z = {1400,1600}({4,6}x,, + 4x15 + 3X13 + 6x14) 
+ {900,1100}(7x5, + 5x29 + 4x93 
+ {3,5} x24 ) 
+ {4500,6000}(4x3, + {9,11 }x3, + 8x33 
+ 5x34 ) 
+ {4000,5000}(6x,, + {2,4}x%42 + 0X43 
+ 0x44 ) + {400,700}(3x., + {10,14}x., 
+ 10x53 + 6x54) ~ Max 


Mathematical model: 
Find the maximum value of 
Z = {1400,1600}({4,6}x,, + 4x15 + 3x13 + 6x44) 
+ {900,1100}(7x51 + 5x22 + 4X23 
+ {3,5} x24 ) 
+ {4500,6000}(4x3, + {9,11 }x3, + 8x33 
+ 5x34) 
+ {4000,5000}(6x,, + {2,4}x%45 + 0X43 
+ 0x44 ) + {400,700}(3x., + {10,14}x5> 
+ 10x53 + 6x54) ~ Max 


Within restrictions: 
X44 + X21 + X34 + X41 + X51 = 60 
X12 + X22 + X32 + X42 + X52 = 150 
X43 + X23 + X33 +X 43 + X53 = 20 


oTATe 
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Xq4 + X24 + X34 + X44 + X54 = 10 
{4,6}x11 + 4x12 + 3x13 + 6x14 = 2500 
7X21 + 5X22 + 4X03 + {3,5}x24, = 1000 
4X3, + {9,11 }x32 + 8x33 + 5x34 = 600 
6X41 + {2,4}x42 + 0X43 + Oxy, > 200 
3x51 + {10,14}xs, + 10x53 + 6x54 > 800 
Xj 20 5i=1,2,3,4,5 and j = 1,2,3,4 
The first issue: 


An executive in one of the companies asked an expert in the 
science of operations research to help him obtain an optimal 
solution through which to achieve the lowest cost of 
transportation and operation of warehouses he wants to 
establish in order to expand the company's work and provided 
him with information through which the expert formulated the 
following issue: The text of the problem according to the 
concepts of neutrosophic science: A retail company plans to 
expand its activities in a specific area by establishing two new 
warehouses, the following table shows the potential locations, 
the number of customers and the possibility of meeting the 
demand for the sites where (+) has been placed in the event that 
the site can meet the customer's request and put (x) the 
opposite and code Nc, For the cost of transferring one unit 


from site t to customer J he got the following table: 
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Customer 
By By B3 By 
site 
Ay * x * 
Ney, Nc42 Ney, 
Az * * * 
Neo Nez NCp3 NCp4 
Az x * * * 
Nc32 Nc33 Nc34 
Customer dD, Dz D3 D4 
orders 


Table (1) Transportation cost in case of location selection 


With the following information available for each of the 
candidate locations for warehouses 


Operating cost Initial Invested Site Capacity 
per unit Capital (Monetary 
(monetary unit) Unit) 
first Np, ky Aj 
second Np2 kp Az 
third Np3 ks A3 


Table (2) operation information 


It is required to choose suitable locations for warehouses that 
make the total costs of investment, operation and transportation 
as small as possible. 


Building the mathematical model: Each site has a fixed capital 
cost independent of the quantity stored in the warehouse 
referred to that site and also has a variable cost proportional to 
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the quantity transported, and therefore the total cost of 
establishing and operating the warehouse is a non-linear 
function of the stored quantity and using binary integer 
variables can be formulated the issue of determining the 
location of the warehouse in a program with integers where we 
assume that the binary integer variabled; Symbolizes the 
decision to choose the site or not to choose it in other words 
pe {t if we chose the site i 
0 otherwise 

Suppose that x,; is the quantity transferred from site & to 
customer J, so the constraint expressing the ability of the first 
site to meet the requests is as follows: 


X11 + X12 + X14 < A,6, 


When 6, = 1, the first location with capacityA, is chosen. The 
quantity transported from the first site cannot exceed the 
capacity of that siteA, When 6, = 0 the non-negative variables 


X41,X12,X14 = 0 directly, indicating that it is not possible to 
ship from the first location 


In a similar way, we obtain the following two constraints for 
the second and third signatories. 


X24 + X22 + X23 + X24 < Az6, 
X34 + X33 + X34 < A363 


To choose exactly two locations, we need the following 
restriction: 


6, +65 + 63 =2 
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As 6, can take one of the values of 0 or 1 only, the new 
constraint will force two variables from among the three 
variables, 6; to be equal to one 


The restrictions for customer requests can be written as 
follows: 


first customer X14 + X21 = D, 

Second customer X12 + X22 + X32 = Dz 
Third customer X23 + x33 = D3 

Forth customer X14 + X24 + X34 = Dy 


To write the objective function, we note that the total cost of 
investment, operation and transportation for the first site is as 
follows: 


Kk, 6, + Npy (X44 + X42 + X44) + NCq1%41 + NCq2%42 
+ NC44X14 


When we do not choose the first site, variable 6, = 0 And 
that forces the variables 


X16 X42° X44 = 0 

In a similar way, the cost functions of the second and third 
sites can be written, and thus the full formulation of the issue 
of assigning the location of the warehouse is reduced to the 
following correct mixed program: Z is meant to be made 
minimal 

Z = ky, + Npy (X44 + X12 + X14) + NCq1%41 + NCq2%X12 + 
NCy4X14 + k252 + Npo(%21 + X22 + X23 + X24) + NCQ1%X21 + 
NC 22X22 + NC23X23 + NC24X24 + k363 + Np3 (X32 + X33 + 
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X34) + NC32X32 + NC33%33 + NC34X34 
considering the following restrictions: 

61 “- b> + 63 as 

X41 +X%21 =D, 

X42 + X22 + X32 = Dz 

X23 + X33 = D3 

6; true variable for i = 1,2,3 
xj 20 ;1=1,2,3 andj = 1,2,3,4 

The second problem: The second request addressed by the 
executive was how I can choose the appropriate projects to 
operate a limited capital available in the company through a 
number of projects presented, through the information provided 


by the official in charge of managing the company, the expert 
formulated the following issue: 


The issue of the capital budget: A company plans to disburse 
its capital during the 7; periods. Where: 


j =1,2,...,n, and there isA;A proposed project where 
it = 1,2,——-—,m versus a limited capital B; Available for 
investment in period j and when choosing any project i 
becomes in need of a certain capital in each period j we 
denote itNaj;. It is a neutrosophic value, the value of each 
project is measured in terms of the liquidity flow 
corresponding to the project in each period minus the value 
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of inflation, and this is called net present value (NPV), we 


denote it Nv; Accordingly, the following table can be 
organized: 
eriod T; T. | weeeennenee T. 
projec 

A, Nay, Nay. ~~ (| ---------- Nay, 

A, Nay, Nay Nagy 

Am NQm1 NQm2 ‘| ---------- Namn 
Limited By By | ---------- Ba 
capital 


Table (3) Return on Investment during Periods 


What is required in this problem is to select the right projects 
that maximize the total value (NPV) of all selected projects. 
Formulation of the mathematical model: 


Here we assume a binary integer variable x,It takes the value 
one if the project J is selected and takes the value zero if the 
project j is not selected 


em 


Then the objective function is given by the following 
relationship: 


if we chose project i 
0 otherwise 
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m 
L.= YY Nv; Xi 
i=1 


Then the objective function is given by the following 
relationship: 


m 
a Najj x; SB; ;sJ=1,..,n 
i=1 


Accordingly, we get the following mathematical model: 


Find the maximum value of the function: 


m 
1. >. Nv; Xi 
i=1 


considering the following restrictions: 


m 

i=1 
x; A binary variable takes one of the values 0 or 1 for all values 
of t = 1,...,m in the previous two issues, we got models with 
integers that have special methods of solution. This research 
cannot be presented and we will present them in later 
research using the concepts of neutrosophic science 


1- Formulation of the problem and the construction of 
mathematical model according to neutrosophic 
values: 


The study concluded in the research [12] shows us how to 
construct neutrosophic linear models, (the linear model is a 
neutrosophic model if at least one of the likes of variables in 
the target function or neutrosophic value constraints) 


The text of the issue: 


2154s 
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The company has n rank for inspectors and wants to assign the 
task of quality control to them, and K pieces should be audited 
daily during an S hour of work per day, in the following table 


we explain the full information about the inspectors and for all 
ranks: 


Number Accuracy | Inspector's Number of | The fine paid 

of pieces | (percent) | remuneration | inspectors | by the 

checked ( Monetary company for 

(hour) Unit per Hour each fault to 

) the inspector 

NM, ND, Gy Ay R 
2 NM, ND, Gz A, R 
n NM, NDy Gn An R 


Table (3 ) Information on inspectors using neutrosophic values 
The number of pieces is a neutrosophic value NM; = M; + &; 
where €; is the indeterminacy on the number of pieces, it can 
take one of the shapes [Aja Ajo or {Ajs, Ajo} or any value close 
to M; as well as the precision, neutrosophic values 
ND; = D, + 6; 
where 6; is the indeterminacy on the precision that can take one 


of the shapes [uj Ljo| or {uj Ljo} or any value close to Dj. 


Required: Formulate the appropriate mathematical model 
through which we can assign the optimal support to the 
inspectors so that the cost of inspection is as low as possible 


Building the neutrosophic mathematical model: 


To build the mathematical model, we impose 
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X1,X2,-— — ——,X, the number of inspectors of each rank on 
the order assigned to the inspection task, then the following 
inequality must be met: 


xj SA; 5 J=H1,2,-——,n 

Since the company needs to audit K piece daily within S 
working hour per day, the following set of restrictions must be 
met: 

n 

>, S(NM;)x; = K 

j=l 
To obtain the target function, we note that the company bears 
two types of costs during the inspection process, the inspector's 
fee and the fine corresponding to the error committed by the 


inspector for each piece then the target follower writes as 
follows: 


100 — ND; 
NZ = G; + (NM; )Rj la Xj 
j=l 
Then the mathematical model is written as follows: 


100 — | 


NZ = s\" G; + (NM,)R; I x; > Min 
j=1 


Within Restrictions 


n 
>, S(NM))x; = K 


j=l 
MeOS fH LAS SHR 
Example 1: 
- 156 - 


Neutrosophic linear models and algorithms to find their optimal solution 


A company has three ranks for inspectors and wants to assign 
the task of quality control to them, and 1500 pieces should be 
audited daily during 8 working hours per day, in the 
following table we explain the full information about the 
inspectors and for all ranks, in this example we will take the 
number of pieces checked by the inspectors from each rank as 


neutrosophic values 


About the | Number of | Accuracy Inspector's Number | The fine paid by 

Inspector pieces (percent) remuneration of the company for 
checked (Monetary Unit | inspectors | each fault to the 

(hour) per Hour ) inspector 
Inspector rank 

1 {15,16} 95 4 10 9 

2 {10,11} 90 3 6 9 

3 {25,26} 98 5 8 9 


Table (4 ) Information on inspectors using neutrosophic values 


Required : Formulate the appropriate mathematical model 
through which we can assign the optimal assignment to the 
inspectors so that the cost of inspection is as low as possible 


To build the mathematical model, we impose x1 x2,X3 as the 
number of inspectors from the three ranks in the order assigned 
to the inspection task, then the following inequality must be 
met: 


Since the company needs to audit K pieces daily within S 
working hour per day, the following set of restrictions must be 
met: 
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n 
> 8M,x; => 1500 
j=l 


That is 
8(M, x, + M2x2 + M3x3) = 1500 
From it we get the following restriction: 
8{15,16}x, + 8{10,11}x, + 8{25,26}x3 = 1500 


To obtain the target function, we note that the company bears 
two types of costs during the inspection process, the inspector's 
fee and the fine corresponding to the error committed by the 
inspector for each piece then the target follower writes as 
follows: 


Then the cost of the inspector is calculated from j the hourly 
rank through the following relationship: 


100 —D, 
) sf=12---m 


From that we get 


100 — 95 
NC, = 4+ {15,16} x 2 x (a) 2765556) 
100 
100 — 90 
NC, = 3 + {10,11} x 2 x (a) = (55.0) 
100 
100 — 98 
NC, = 5 + {25,26} x 2 x (—-} = (6,6.04} 


The total costs for all inspectors assigned to the task of quality 
control per hour shall be given by the following relationship: 


NTC = [ +(NM;)R (=) 
p> gessa ek ® OME ago j 
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NTC, = {5.5,5.6}x, + {5,5.2}x2 + {6,6.04}x; 


substituting the following target phrase: 


NZ >. [ +(NM,)R —] 
= . . . EO X; 
24 J JZ°J 100 J 


We get: 
NZ = {44,44.8}x, + {40,41,6}x5 + {48,48.32}x3 
From the above, we can develop the following mathematical 
model: 
We want to find : 
Min(NZ) = {44,44.8}x, + {40,41,6}x, + {48,48.32}x, 
Within Restrictions 
x, < 10 
Xz < 6 


x3 <= 8 


8{15,16}x, + 8{10,11}x, + 8{25,26}x, > 1500 
x, 20; fj=1,2,3 
Example 2: 


A company has three ranks for inspectors and wants to assign 
the task of quality control to them, and 1500 pieces should be 
checked daily during 8working hours per day, in the following 
table we explain the full information about inspectors and for 
all ranks, in this example we will take the accuracy of 
inspection for each inspector as neutrosophic values in the 
form of areas whose minimum range is less accurate and the 
highest range is the highest accuracy that the inspector reaches 
by rank 
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About the Number of | Accuracy Inspector's Number : : 
: ; The fine paid by 
pieces (percent) remuneration of 
hecked ; ; the company for 
checke (monetary unit inspectors sac fault tothe 
(hour) per hour) ? 
inspector 
Inspector rank 
1 15 [95,97] 4 10 2 
2 10 [90,92] 3 6 9 
3 25 [98,99.5] 5 8 9 


Table (5) Information on inspectors using neutrosophic values 


Required : Formulate the appropriate mathematical model 
through which we can assign the optimal assignment to the 
inspectors so that the cost of inspection is as low as possible 


To build the mathematical model, we impose xx 2,X3 the 


number of inspectors from the three ranks in the order assigned 
to the inspection task, then the following inequality must be 


met: 


IA IA 


IA 


10 
6 
8 


Since the company needs to audit K pieces daily within S 
working hour per day, the following set of restrictions must be 


met: 


That is 


n 
> 8M;x; => 1500 
j= 


8(M, x, + M2x2 + M3x3) = 1500 
- 160 - 


Neutrosophic linear models and algorithms to find their optimal solution 


From it we get the following entry: 
120x, + 80x2 + 200x3 = 1500 


To obtain the target function, we note that the company bears 
two types of costs during the inspection process, the inspector's 
fee and the fine corresponding to the error committed by the 
inspector for each piece then the target follower writes as 
follows: 


Then the cost of the inspector is calculated from j the hourly 
rank through the following relationship: 


100 — ND; 
NG = Gj + M)Rj)(——) ; j=1,2,-——--,n 
From that we get 
NC ae ge (eee [4.9,5.5] 
oo 100 atria 
100 — [90,92] 
NG, = 3+ 10% 2%|-—— = [4.6,5] 


100 — [98,99.5] 


NC3=5+25x2x 
are ( 100 


= [5.25,6] 


The total costs for all inspectors assigned to the task of quality 
control per hour shall be given by the following relationship: 


NTC > " [ +M;R (= 2 — 
ati, A é ‘e ——— X; 
J iA J Jed 100 J 


NTC, = [4.9,5.5]x, + [4.6,5]x2 + [5.25,6]x; 


substituting the following target phrase: 


NZ = a) rae | fa (enc 
{4 js ap9 100 ‘| 
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We get: 
NZ = [39.2,44]|x, + [36.8,40]x, + [42,48]x, 
From the above, we can develop the following mathematical 
model: 
We want to find : 
Min(NZ) = [39.2,44]x, + [36.8,40]x, + [42,48]x, 

Within Restrictions 

x, < 10 

Xz < 6 

x35 8 

120x, + 80x2 + 200x3 = 1500 
xj 20; fj=1,2,3 

In the two examples! and 2 for the optimal solution we use the 
neutrosophic simplex method sufficiently explained in the 
reference [13]. 
From the previous model, we notice that x; takes a positive 
value only when oF = 1, and in this case, the production of the 
product j is limited by the quantity d; and the fixed production 
cost K; is included in the goal function 
The idea of indeterminacy is the basis of neutrosophic science 
represented here through the use of the binary integer variable 
because the optimal solution depends on the decision to 
produce a product or not to produce it, although we cannot 
guarantee the company a safe working environment in light of 


the great changes in the labor market through price strikes, 
resource availability or non-availability, and so on. 
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So it was necessary to reformulate this problem using 
neutrosophic values for the sales opportunity d; and the cost of 
producing one unit of each product Cjand selling price P; so 
that the sales opportunity becomes d; + € , production cost 
C; +; and selling price Pj) + p; where €; andy; and @; are 
the indeterminacy and is the change in the sales opportunity, 
cost and selling price respectively depending on the conditions 
of the work environment and takes one of the following forms: 


Which are values close to the values d; and C; and can be any 
neighborhood to them. 

Then the text of the problem becomes as follows: 

The text of the problem according to neutrosophic science: 
A company is planning to produce N product where the 
product j! needs a fixed preparation cost or a fixed production 
cost Kj; independent of the quantity produced, and needs a 
variable cost C; + 4; per production unit commensurate with 
the quantity produced, we suppose that each unit of the product 
j needs a;; a unit of the supplier t where there is M supplier. 
Assuming that the product j that has a sales opportunity d; + & 
is sold at the price of P; + @; monetary unit per unit and that 
only b; unit of the supplier i is available where i = 1,2,—— 
—,M_ the goal of the problem becomes to determine the 


optimal product mix that makes the net profit as great as 
possible. 
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Formulation of the mathematical model: 
Determination of the cost: 


From the text of the problem, we note that the total cost of 
production is the fixed cost in addition to the variable cost, 
which is a nonlinear function of the quantity produced, 

But with the help of binary integer variables 6,, the problem 
can be formulated in the form of a linear model with integers 
We assume that the binary integer variable 6; symbolizes the 
decision to produce the product j or not to produce it in other 


words 


pte ( j if production desicion was taken 
a 0 otherwise 


Then the cost of producing one unit of the product becomes as 
follows Kj;6; + (C; + uj;)x;, where 6; = 1 if x; > 0 and 6; = 0 
if x; = 0 and therefore the goal function becomes as follows: 


N N 
aS yy pia) i » UGG + G + Hx) 
j= iz 
Restrictions of the problem: 


A restriction on the supplier i is given in the following 
relationship: 


N 
> Ajj Xj Sd; ‘tL=1,2,-—--—,M 
j= 
The restriction of the demand for the product j is given by the 

following relationship: 


xj < (d; + Ej) 6; J = 1,2,- = —,N 
Mathematical model: Find the maximum value of the function: 
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N N 
Z= >». Oy -\ (G5 + (Cj + H,)x;) 
j= j= 


Within Restrictions 


N 

>. Ajj Xj Sd; ‘i=1,2,-——--,M 
j=t 

xj < (dj + €)6; j= 1,2,-—-—-—,N 


x; 2 0 and 6; = 1 or 6; = 0 
And or for all values 

J = 152; ee ee (ad N 
From the previous model, we note that x; takes a positive value 
only when 6; = 1 and in this case the production of the product 
j is limited by the quantity d; + €; and the fixed production 
cost K; included in the goal function, by solving this model 
we get an optimal neutrosophic value for the goal function 
NZ* through which we know the profit that the company can 


achieve in the best and worst conditions and enable the 
company to develop appropriate plans for the workflow in it. 
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Conclusion: 


Through the previous study, we note that by using the linear 
programming method, we can provide the optimal solution to 
most of the problems that production companies can face by 
formulating the situation under treatment with a problem that 
can be converted into a linear model, the optimal solution of 
which helps decision-makers make optimal decisions for the 
workflow so that the greatest profit is achieved. It is possible, 
and in order to obtain solutions that have a margin of freedom, 
we can use the data, neutrosophic values, values that take into 
account all the conditions that the system represented by the 
linear model may experience. 
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